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ABSTRACT 


Title  of  Dissertation:  THE  EFFECTS  OF  NOISE  ON 
NONLINEAR  SYSTEMS  NEAR  CRISIS 
John  Curtis  Sommerer,  Doctor  of  Philosophy,  1991 

Dissertation  directed  by  Professor  Edward  Ott,  Ph.  D., 

Department  of  Physics  and  Astronomy,  and 
Department  of  Electrical  Engineering 

We  consider  the  influence  of  random  noise  on  low-dimensional,  nonlinear 
dynamical  systems  with  parameters  near  values  leading  to  a  crisis  in  the  absence  of  noise. 
In  a  crisis,  one  of  several  characteristic  changes  in  a  chaotic  attractor  takes  place  as  a 
system  parameter  p  passes  through  its  crisis  value  pc.  For  each  type  of  change,  there  is  a 
characteristic  temporal  behavior  of  orbits  after  the  crisis  (p>pc  by  convention),  with  a 
characteristic  time  scale  t.  For  an  important  class  of  deterministic  systems,  the 
dependence  of  Ton p  is  t ~(p-pc)~Y for  p  slightly  greater  than pc-  When  noise  is  added  to  a 
system  with  p<pc,  orbits  can  exhibit  the  same  sorts  of  characteristic  temporal  behavior  as 
in  the  deterministic  case  for  p>pc  (a  noise-induced  crisis).  Our  main  result  is  that  for 
systems  whose  characteristic  times  scale  as  above  in  the  zero-noise  limit,  the 
characteristic  time  in  the  noisy  case  scales  as  t~a~Y  g[(pc  -p)/o],  where  <x  is  the 
characteristic  strength  of  the  noise,  g(«)  is  a  non-universal  function  depending  on  the 
system  and  noise,  and  yis  the  critical  exponent  of  the  corresponding  deterministic  crisis. 
Illustrative  numerical  examples  are  given  for  two-dimensional  maps  and  a  three- 
dimensional  flow.  In  addition,  the  relevance  of  the  noise  scaling  law  to  experimental 
situations  is  discussed. 


We  investigate  experimentally  the  scaling  of  the  average  time  f  between 
intermittent,  noise-induced  bursts  for  a  chaotic  mechanical  system  near  a  crisis.  The 
system  studied  is  a  periodically  driven,  variable-noise,  magnetoelastic  ribbon.  We 
determine  /for  the  low-noise  (“deterministic”)  system,  then  add  noise  and  observe  that 


the  scaling  for  ris  as  predicted. 
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I.  INTRODUCTION 


Crises1  are  commonly  seen  in  deterministic,  dissipative,  nonlinear 
dynamical  systems,  and  they  have  been  demonstrated  in  many  experimental2,3  and 
numerical  studies.  In  a  crisis,  one  of  several  characteristic  changes  in  a  chaotic  attractor 
takes  place  as  a  system  parameter  p  passes  through  some  critical  value  pc,  at  which  the 
crisis  occurs:  the  chaotic  attractor  may  be  destroyed,  it  may  increase  in  size,  or  it  may 
merge  with  some  other  chaotic  attractor.  As  detailed  in  the  next  paragraph,  typical 
trajectories  of  a  post-critical  system  (by  convention  p  >  pc)  behave  like  those  of  the  pre- 
critical  system  for  a  time  that  is  sensitively  dependent  on  initial  conditions.  For  p>Pc  and 
randomly  chosen  initial  conditions  in  the  basin  of  the  pre-critical  p<Pc  attractor,  the 
length  of  this  time  T  is  found  to  be  exponentially  distributed,4  with  a  mean  length  x  that 
depends  on  p;  i.e.,  the  probability  density  of  T  is  PT{t)-X~^ertl‘T,  t> 0. 

The  definition  of  the  characteristic  time  x  differs  for  the  three  types  of 
crises  mentioned  above.  In  the  case  of  attractor  destruction,  the  trajectory  initially  moves 
chaotically  as  though  it  were  still  confined  to  the  pre-critical  attractor,  and  then  leaves  the 
region  in  phase  space  where  the  former  attractor  was  located,  and  tends  to  some  other 
attractor,  or  possibly  to  infinity.  The  period  of  time  that  the  post-critical  system  behaves 
like  the  pre-critical  system  is  referred  to  as  a  chaotic  transient.  For  this  case  x  is  the  length 
of  the  transient,  averaged  over  many  different  initial  conditions.  In  the  case  of  attractor 
enlargement,  the  post-critical  trajectory  stays  in  the  phase  space  region  of  the  the  pre- 
critical  chaotic  attractor  for  some  time,  and  then  bursts  into  chaotic  motion  over  a  larger 
region  in  phase  space,  before  being  reinjected  into  the  region  of  the  pre-critical  attractor, 
the  process  then  repeats.  For  this  case,  x  is  the  average  time  between  such  bursts.  In  the 
case  of  attractor  merging,  which  can  occur  in  systems  with  symmetries,  the  post-critical 
motion  switches  repeatedly  between  the  phase-space  regions  occupied  by  the  pre-critical 
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chaotic  attractors  that  merge  at  the  crisis.  The  characteristic  time  T  is  now  the  average 
time  between  switchings. 

For  a  large  class  of  low-dimensional  systems  (including  one-dimensional 
maps,  and  many  two-dimensional  maps  and  three-dimensional  flows),  the  characteristic 
time  T  is  found  to  scale  with  the  parameter  as5 

t~\p-pc\-y.  (i) 

The  crises  with  such  scaling  behavior  result  from  the  collision  of  the  chaotic  attractor 

with  the  stable  manifold  of  an  unstable  periodic  orbit.  In  particular,  when  we  consider 

two-dimensional  maps  or  three-dimensional  flows,  we  are  interested  in  cases  where  the 

collision  is  a  tangency  between  the  stable  manifold  of  an  unstable  periodic  orbit  and  the 

unstable  manifold  of  the  same  (homoclinic  tangency  crisis)  or  another  (heteroclinic 

tangency  crisis)  unstable  periodic  orbit.  The  critical  exponent  y  in  the  scaling  law  (1)  is 

determined  by  the  eigenvalues  of  the  periodic  orbit  whose  unstable  manifold  is  involved 
in  the  tangency^  (or  is  ^  for  one-dimensional  maps  with  generic,  quadratic  maxima). 

In  a  deterministic  system  for  p<Pc>  the  characteristic  time  is  infinite, 
because  the  orbit  remains  forever  on  the  pre-critical  attractor.  However,  if  some  random 
noise  is  added  to  the  system,  there  is  the  possibility  that  a  trajectory  starting  in  the  basin 
of  attraction  of  the  deterministic  system’s  chaotic  attractor  will  behave  like  an  orbit  of  the 
post-critical  system  even  for  p<Pc ,  thereby  producing  a  transient  response.  For  example, 
in  the  cases  of  attractor  destruction  and  merging,  noise  can  kick  an  orbit  across  the 
boundary  of  the  attractor’s  basin  of  attraction.  This  situation  is  referred  to  as  a  noise- 
induced  crisis,  and  it  stands  in  contrast  to  those  metastable  situations  where  the 
corresponding  deterministic  system  is  neither  chaotic  nor  near  crisis.  Because  real 
physical  systems  are  always  accompanied  by  noise,  we  believe  that  noise-induced  crises 
may  be  an  important  complication  in  the  experimental  investigation  of  transient  chaos. 
Alternatively,  the  intentional  addition  of  noise  to  an  experiment  may  offer  another 
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dimension  of  control  in  understanding  the  underlying  crisis  dynamics.  It  is  the  purpose  of 
this  dissertation  to  provide  a  theory,  supported  by  numerical  and  physical  experiments, 
for  the  effects  of  noise  on  crises. 

The  production  of  metastable  states  in  dynamical  systems  by  the 
introduction  of  noise  has  been  considered  in  the  literature  of  physics  and  stochastic 
processes  for  many  years.6  The  usual  approach  has  been  to  describe  the  evolution  of  a 
probability  density  by  a  partial  differential  equation,  obtaining  state  lifetimes  by 
determining  the  rate  at  which  probability  mass  passes  through  some  boundary  in  phase 
space.  Such  an  approach  is  in  principle  applicable  to  noise-induced  crises,  but  the 
underlying  nonlinear  dynamics  typically  make  for  an  intractable  density-evolution 
equation. 

The  effects  of  noise  on  transitions  to  chaos  via  the  period-doubling  and 
intermittency  routes  have  been  considered  previously.7  Numerical  studies  of  noise- 
induced  crises  have  also  been  performed. In  Ref.  9,  which  included  numerical 
experiments  and  an  analysis  of  a  one-dimensional  map,  it  was  suggested  that  noise- 
induced  crises  may  admit  scaling  behavior  for  their  characteristic  times.  In  this 
dissertation,  we  argue  that  for  noise-induced  crises,  corresponding  to  deterministic  crises 


with  scaling  behavior  (1),  the  characteristic  time  scales  as 

t~  a-y  (2) 

where  <r  is  a  measure  of  the  strength  of  the  noise,  g(*)  is  a  non-universal  function 
depending  on  the  system  and  the  distribution  function  of  the  noise,  and  y  is  the  critical 
exponent  of  the  corresponding  deterministic  crisis.  We  obtain  this  result  without  recourse 
to  density-evolution  partial  differential  equations.  [Note  that  consistency  of  (1)  and  (2) 
implies  that  g(-x)~x  -f  for  large  positive  x.] 

In  Chapter  II,  we  derive  scaling  law  (2),  starting  with  one-dimensional 
maps  and  generalizing  to  higher  dimensions.  We  also  review  the  deterministic  theory 
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leading  to  scaling  law  (1),  because  it  determines  the  critical  exponent  in  the  new  scaling 
law  (2),  and  underlies  certain  features  seen  in  the  numerical  examples  with  noise 
discussed  later. 

Chapter  III  presents  a  more  detailed  analysis  of  the  one-dimensional  map 
discussed  in  Ref.  9.  In  determining  the  function  #(•)  of  Eq.  (2)  explicitly  for  this  case 
(something  we  are  unable  to  do  for  the  more  complicated  examples),  we  emphasize  how 
certain  simplifying  assumptions  made  in  the  derivation  of  Eq.  (2)  affect  the  accuracy  and 
applicability  of  the  scaling  law.  We  then  outline  the  procedures  used  to  get  the  numerical 
results  in  the  balance  of  this  dissertation.  Numerical  results  for  the  one-dimensional  map 
show  that  Eq.  (2)  is  satisfied  for  this  case  (see  also  Ref.  9),  and  support  the  exponential 
distribution  of  t. 

Chapter  IV  includes  two  examples  of  two-dimensional  maps,  which  have 
deterministic  crises  with  critical  exponents  depending  on  the  eigenvalues  of  mediating 
unstable  saddles.  One,  with  a  heteroclinic  tangency  crisis,  is  used  to  demonstrate  the 
effect  of  the  noise  distribution  function  on  the  function  #(•)•  The  other  system 
demonstrates  that  (2)  works  equally  well  for  a  homoclinic  tangency  crisis. 

In  Chapter  V,  we  discuss  two  potential  complications  to  observing  scaling 
law  (2).  In  the  first  example,  we  discuss  circumstances  that  produce  results  clustering 
around,  rather  than  on,  the  curve  given  by  the  function  #(•).  A  second  example  shows 
how  the  presence  of  noise  can  lead  to  trajectories  escaping  via  routes  different  from  those 
of  the  corresponding  deterministic  crisis;  such  alternate  routes  of  escape  produce  a  more 
stringent  bound  on  a  than  would  otherwise  be  required  for  scaling  to  be  observed. 

Although  the  derivation  of  scaling  law  (2)  in  Chapter  II  considers  only 
maps,  in  Chapter  VI  we  discuss  its  applicability  to  flows  as  well,  and  present  a  supporting 
numerical  example  in  which  the  noise  is  introduced  in  a  physically  motivated  way  to 
produce  a  stochastic  differential  equation. 
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In  Chapter  VII,  we  discuss  the  implications  of  Eq.  (2)  for  experiments 
studying  transient  chaos,  and  in  Chapter  VIII,  we  describe  a  physical  experiment 
confirming  for  the  first  time  the  applicability  of  Eq.  (1)  [specifically,  the  prediction  of  y 
in  Eq.  (1)]  and  Eq.  (2)  to  physical  systems.  The  final  chapter  briefly  summarizes  our 
conclusions.  Three  appendixes  contain  technical  details  of  the  derivation  of  g(*)  for  the 
logistic  map,  of  the  integration  schemes  used  to  study  the  stochastic  differential  equation, 
and  of  the  estimation  of  eigenvalues  of  unstable  periodic  orbits  from  experimental  time 
series. 
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H.  DERIVATION  OF  SCALING  LAW  (2) 


A.  One-dimensional  maps 

As  discussed  in  Chapter  I,  we  are  interested  in  crises  which  result  from  the 
collision  of  a  chaotic  attractor  with  the  stable  manifold  of  an  unstable  periodic  orbit.  For 
concreteness,  consider  a  crisis  in  one  dimension  mediated  by  an  unstable  fixed  point.  The 
system  is  governed  by  the  one-dimensional  map 

xn+i  =F(xn;p)  +  Zn,  (3) 

where  F(»)  is  the  deterministic  map,  and  the  Zn  are  independent,  identically  distributed 
random  variables  with  density  function  P^z).  The  deterministic  crisis  occurs  when,  as  the 

parameter  is  increased  through  pc,  the  chaotic  attractor  (an  interval)  collides  with  the 
unstable  fixed  point.  We  treat  only  nondegenerate  crises,  that  is  those  where  the  distance 
e  from  the  attractor  to  the  mediating  unstable  fixed  point  scales  linearly  with  Pc~P  for  p 
near  pc.  The  noise-induced  crisis  results  from  noise  kicking  a  trajectory  across  the  e-gap 
(or  its  preimage)  between  the  attractor  and  the  mediating  unstable  fixed  point  (or  its 
preimage). 

We  estimate  the  rate  at  which  such  events  occur,  using  several  simplifying 
assumptions  to  apply  random  variable  techniques.  The  procedure  is  illustrated 
schematically  in  Fig.  1.  First,  we  suppose  that  the  invariant  density  for  the  deterministic 
attractor  at  p  =pc  does  not  change  in  form  as  p  is  decreased  slightly,  but  rather  that  the 
attractor  measure  simply  shifts  away  from  the  fixed  point  linearly  with  pc-p •  In  Fig.  1(a) 
the  attractor  measure  is  schematically  represented  as  a  smooth  function  which  is  zero 
within  a  distance  e~  pc  —  p  from  the  fixed  point.  We  further  assume  that  iterates  of  the 
map  can  be  treated  as  independent  realizations  of  a  random  variable  with  density  function 
equal  to  the  invariant  density  of  the  deterministic  attractor.  Obviously,  this  assumption 
ignores  the  dynamics  near  the  fixed  point.  However,  the  analysis  of  the  logistic  map 
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given  in  the  next  section  will  show  how  consideration  of  the  detailed  dynamics  in  the 
immediate  vicinity  of  the  fixed  point  leaves  the  form  of  the  scaling  law  unchanged.  We 
therefore  only  require  sequences  of  iterates  in  the  vicinity  of  the  fixed  point  (“close 
approaches”)  to  be  sufficiently  separated  in  time  and  sufficiently  uncorrelated  to  provide 
a  Poisson  process  for  the  close  approaches.  But  this  is  only  the  same  requirement  as  for 
scaling  law  (1)  to  apply,  and  is  satisfied  by  the  chaotic  dynamics. 


unstable 
fixed  point 


attractor  measure 
in  local  coordinates  \ 


finite-variance 
/  noise  density 


(a) 


convolution  (b) 

I 


unstable ,  kPx+zM 
fixed  point 


fraction  of 
trajectories 
escaping 


- £  0 

dynamical  density 

(c) 


FIG.  1.  The  natural  measure  of  the  attractor,  here  drawn  as  smooth,  is  treated  like  the 
density  of  a  random  variable  and  convolved  with  the  noise  density  to  produce  the 
dynamical  density  from  which  the  escape  rate  is  calculated. 
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We  account  for  the  addition  of  noise  to  the  system  by  convolving  the 
invariant  density  of  the  attractor  with  the  density  function  of  the  noise  [Fig.  1(b)]  to 
construct  a  “dynamical  density”  [Fig.  1(c)].  Note  the  implicit  assumption  that  the  additive 
noises  ...,  Zrt.3,  Zn.2,  Zn.i  have  no  effect  on  the  density  of  xn.  We  associate  the 
probability  mass  falling  on  the  other  side  of  the  unstable  fixed  point  from  the 
deterministic  attractor  with  those  trajectories  that  terminate  the  transient  response  of  the 
metastable  system.  Thus  we  calculate  the  characteristic  time  for  the  noise-induced  crisis 
as  proportional  to  the  reciprocal  of  the  escaping  probability  mass.  This  ignores  orbits  that 
escape  from  the  attractor  and  are  then  reinjected  by  the  noise  (as  well  as  all  higher-order 
contributions),  but  as  we  will  show  in  Chapter  III,  the  linear  scaling  of  the  close  approach 
region  on  forward  or  backward  iteration  of  the  map  means  that  all  such  contributions  will 
scale  the  same  way,  and  hence  do  not  alter  the  form  of  Eq.  (2). 

We  make  two  further  assumptions  to  derive  scaling  form  (2)  via 
dimensional  analysis.  First,  we  assume  that  the  invariant  density  of  the  attractor, 
expressed  in  the  local  coordinates  of  Fig.  1  and  for  small,  positive  x  can  be  written  as 

p( x)  -  x?~l  [for  *<0,  p(jc)=0)]  For  one-dimensional  maps,  the  critical  points  [i.e.,  the 
mavima  and  minima  of  F(*)]  are  genetically  quadratic.  This  leads  to  y  being  genetically  2 

for  one-dimensional  maps.  For  non-generic  cases  when  x  is  near  a  critical  point  xc,  we 
have  F(x)~  F(xc)+K\x-xc\s,  (5*2),  and  the  exponent  y  is10  1/5.  We  therefore  leave  y 
general.  Second,  we  assume  that  the  density  of  the  additive  noise  P ^z)  falls  off  rapidly 

enough  with  z  that  various  limit  and  asymptotic  substitutions  in  the  following  derivation 

introduce  negligible  error.  The  reciprocal  of  the  characteristic  time  f  is  proportional  to  the 

shaded  probability  mass  in  Fig.  1(c),  which  is  the  definite  integral  of  a  convolution: 

—£  —£  00 

1/t  ~  J  dw  Fx+Z(w)  =  |  dw  jdx  p(x)  Pz(w-x). 

_ OO  —OO  —OO 
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Substituting  the  asymptotic  form  of  p(x )  and  writing  Pz(z)  in  the  natural  form 
<t—1  q(z/o),  where  <7  is  the  scale  length  of  the  noise,  we  have: 

1/t  ~  \dw  jdx  ^q €Kx)xy~l , 

— oo  — oo  V  / 


(where  &(•)  denotes  the  unit  step  function).  Making  the  change  of  variables  y=x/<7,  and 

doing  the  second  integral  by  parts  (the  surface  terms  vanish  because  the  noise  density 

vanishes  at  infinity  and  the  invariant  density  in  Fig.  1(a)  vanishes  for;t<0),  we  have: 

-e  oo 

1/t  ~  cV-1  [  dw  [  dy  q'(w/o-y)  yY. 

-oo  0 


Naming  the  second  integral  (a  function  only  of  w/cf)  r(*)  and  making  the  final  substitution 
s=w/a  yields: 

-eJo 


Jtfcr(r). 

— OO 


Naming  the  remaining  integral  (a  function  only  of  eld)  1  /#(•)>  we  obtain  scaling  form  (2). 
Thus  the  critical  exponent  of  the  deterministic  system,  fixed  by  the  properties  of  the 
measure  of  the  near-critical  attractor,  also  governs  the  transient  response  of  the  noisy 
system. 
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B.  Two-dimensional  maps 


Noisy  two-dimensional  maps  of  the  form 

xn+i  =F(x„;p)  +  Zn,  (30 

where  now  x,  F,  and  Z  are  vector-valued,  can  also  be  considered.  This  situation  is  more 
complicated,  because  at  crisis  the  unstable  manifold  of  the  mediating  fixed  point  has  an 
infinite  number  of  tangencies  with  the  stable  manifold  with  which  it  collides. 

Furthermore,  these  tangencies  may  have  different  orientations  in  phase  space,  and  the 
noise  density  Fz(z)  may  be  anisotropic.  However,  the  various  tangencies  are  again  locally 

only  linear  rescalings  of  one  another.  By  projecting  both  the  attractor’s  invariant  measure 
and  the  noise  density  onto  an  axis  normal  to  the  stable  manifold  in  the  vicinity  of  each 
tangency  (as  illustrated  in  Fig.  2),  the  analysis  for  one-dimensional  systems  can  be 
applied  as  before,  where  the  local,  one-dimensional  coordinate  of  the  stable  manifold 
takes  the  place  of  the  unstable  fixed  point  in  the  one-dimensional  case.  Noise  anisotropy 
may  change  the  contributions  of  the  various  tangencies  to  the  total  escape  rate,  relative  to 
isotropic  noise,  but  should  not  typically  change  scaling  form  (2),  where  now  y  is  fixed  by 
the  properties  of  the  measure  of  the  higher-dimensional  chaotic  attractor. 

A  difference  from  the  one-dimensional  situation  is  that  the  critical 
exponent  y  is  no  longer  genetically  confined  to  be  Reference  5  derived  formulas  for 

determining  the  critical  exponent  governing  a  deterministic  crisis  in  a  two-dimensional 
map.  These  derivations  are  reviewed  here,  because  they  help  to  explain  several  aspects  of 
the  numerical  examples  to  be  discussed  in  later  sections. 
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\ 

\ 


line  tangent  to  attractor  A 
and  parallel  to  stable  manifold 


FIG.  2.  For  two  dimensional  maps,  the  natural  measure  of  the  attractor  is  projected 
into  one  dimension  running  normal  to  the  mediating  stable  manifold,  before  being 
convolved  with  the  similarly-projected  noise  density.  Region  R  is  a  rectangular  strip 
of  width  x,  with  only  those  parts  of  its  boundary  that  are  parallel  to  the  stable 
manifold  intersecting  the  attractor. 
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The  tangency  causing  the  crisis  can  occur  in  two  possible  ways. 

(i)  Heteroclinic  tangency.  In  this  case,  the  stable  manifold  of  an  unstable 
periodic  orbit  (B)  is  tangent  to  the  unstable  manifold  of  an  unstable  periodic  orbit  (A)  on 
the  attractor,  as  in  Fig.  3(a).  (Figure  3  is  drawn  for  the  case  of  fixed  points;  for  the  case 
where  A  and  B  are  periodic  points,  similar  pictures  can  be  drawn.) 

(ii)  Homoclinic  tangency.  In  this  case,  the  stable  and  unstable  manifolds 
of  an  unstable  periodic  orbit  (B)  are  tangent,  as  in  Fig.  3(b). 


FIG.  3.  (a)  Schematic  illustration  of  heteroclinic  tangencies  of  the  stable  manifold  of 
the  unstable  periodic  orbit  B  and  the  unstable  manifold  of  the  unstable  periodic  orbit 
A.  (b)  Schematic  illustration  of  homoclinic  tangencies  of  the  stable  and  unstable 
manifolds  of  the  unstable  periodic  orbit  B.  The  figure  is  drawn  for  period-1  orbits. 
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As  in  the  one-dimensional  case,  we  assume  that  the  crisis  is 
nondegenerate,  and  further  that  the  tangency  is  quadratic.  In  both  cases,  at  p~Pc  the 
chaotic  attractor  is  the  closure  of  one  of  the  branches  of  the  unstable  manifold  of  B.  In  the 
heteroclinic  case,  the  critical  exponent  is  given  by 

y=|+/nlail(l/nla2H)-1,  (4) 

where  a\  and  a2  are  the  expanding  (I «i I>1)  and  contracting  (la2l<l)  eigenvalues, 
respectively,  of  the  periodic  orbit  A  in  Fig.  3(a).  In  the  case  of  a  homoclinic  tangency,  we 
have 

y=  lrt\pi\  /  ln\P\Pi\2,  (5) 

where  Pi  and  pi  are  the  expanding  and  contracting  eigenvalues  of  the  periodic  orbit  B  in 
Fig.  3(b).  For  both  cases,  letting  the  contraction  rate  be  large  (la2l-1,  l/fel-1  °°)  recovers 
the  generic  exponent  for  one-dimensional  maps,  7=5. 

For  the  heteroclinic  crisis,  as  p  is  increased  past  pc,  the  unstable  manifold 

of  A  crosses  the  stable  manifold  of  B  (cf.  Fig.  4).  Before  the  crisis  the  attractor  was 

confined  to  the  region  to  the  right  of  the  upper  stable  manifold  segment  of  B,  and  B  is  not 

on  the  attractor  After  the  crisis  an  orbit  initially  in  the  region  in  which  the  chaotic 

attractor  was  confined  for  p<pc  can  eventually  land  in  the  shaded  region  ab  of  Fig.  4. 

Such  an  orbit  will  then  be  attracted  along  the  stable  manifold  of  B  and  then  rapidly  leave 

the  region  to  which  the  chaotic  attractor  was  confined  along  the  left  branch  of  the 

unstable  manifold  of  B.  For  p  near  pc,  the  dimensions  of  region  ab  are  of  order  r  and  r^2, 

where  r~p—pc.  We  now  iterate  the  region  ab  backwards  in  time  for  n  steps.  For  small  r 

and  large  enough  (but  not  too  large)  n,  except  for  the  first  few  backward  iterates,  the 

change  in  the  region  ab  is  governed  by  the  linearization  of  the  map  about  A.  Thus  the 
preiterated  region  a'b'  has  dimensions  of  the  order  r/a^  and  r^/a”,  as  shown  in  Fig.  4. 
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FIG.  4.  Schematic  diagram  illustrating  the  derivation  of  Eq.  (4) 


Because  after  falling  in  region  a'b'  the  orbit  soon  falls  in  region  ab,  we  estimate  T  as  the 

average  time  it  takes  an  orbit  to  land  in  region  a'b' .  Now  consider  the  probability  measure 

of  the  attractor  at  p=Pc •  The  quantity  t  is  then  estimated  as  the  probability  that,  on  a 

given  iterate,  an  orbit  on  the  p=Pc  attractor  falls  in  the  region  a'b',  and  we  denote  this 

probability  by  p(r).  Now  reduce  r  by  the  factor  ai  and  consider  the  resulting  region  ab. 

After  we  iterate  backwards  n+ 1  steps  (instead  of  ri),  the  long  dimension  of  the  preiterated 
region  is  again  r!a\  but  the  width  is  changed  to  (a2^)1/,2/«”+1.  Assuming  that,  for  our 


purposes,  the  attractor  measure  can  be  treated  as  if  it  were  smooth  in  the  direction  of  the 


unstable  manifold  of  A,  we  have 


Mr) 

P(cc2r) 


rW/ax 

(a2r)ll'1laxl 


a\ 


With  the  assumption  that  p(r)~ry,  Eq.  (4)  then  follows.  [Actually  p(r)  is  typically  not 
smooth.  For  the  purposes  of  this  dissertation  it  suffices  to  treat  it  as  if  it  is  smooth.] 


For  the  homoclinic  case,  consider  the  situation  at  p=Pc  represented 
schematically  in  Fig.  5.  We  denote  the  measure  of  the  attractor  in  the  shaded  region 
defined  by  the  unstable  manifold  segment  aoc  and  the  vertical  line  abc  by  p{e),  where  we 
take  the  vertical  line  abc  to  be  distance  £  from  the  stable  manifold  of  B.  We  assume  that 
p.(e)~£Y,  for  small  e,  and,  as  above,  identify  the  exponent  for  p(£)  with  that  governing  the 
scaling  of  the  characteristic  time.  The  basis  for  this  assumption  is  that  for  p  slightly 
greater  than  pc,  we  expect  the  unstable  manifold  of  B  shown  to  poke  over  to  the  other  side 
of  the  stable  manifold  by  a  distance  £  ~  p-pc,  as  in  Fig.  4.  Thus  we  have 


FIG  5.  Schematic  diagram  illustrating  the  derivation  of  Eq.  (5) 
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Also  shown  in  Fig.  5  is  another  vertical  line  segment  def  which  has  been 
chosen  a  distance  ejh.  from  the  stable  manifold  of  B.  The  measure  of  the  attractor  in  the 
region  defod  is  fiiefii)-  Our  goal  in  what  follows  will  be  to  estimate  the  ratio 
By  doing  this  we  shall  be  able  to  determine  the  exponent  y. 

Imagine  that  we  iterate  the  shaded  region  backward  in  time  many  iterates 
so  that  the  points  (a,  b,  c,  d,  e,f,  o)  map  to  ( a',  b',  c',  <£,  e',f ,  o')  which  are  close  to  B  (for 
e  small).  Now  iterate  the  region  d'e'f o' d'  backward  one  further  iterate  to  <f'e"f'o"d" . 
Because  the  primed  letters  are  close  to  B  this  one  further  backward  iterate  is  governed  by 

the  linearized  map  at  B  (i.e.,  by  the  eigenvalues  /?i  and  f$2  evaluated  at  B).  Thus  the 
distance  from  o'  to  e  is  stretched  by  (3  2  t0  o"e",  while  segment/ o' d'  is  compressed  by 

to  become  the  shorter  segment  f'o"d" .  Note,  however,  that  because  the  original 

distance  from  o  to  e  (namely,  £^2)  was  chosen  to  be  shorter  by  precisely  the  factor  P2 
than  the  distance  from  otob  (namely,  e),  we  now  have  that  the  distance  from  e"  to  o"  is 
the  same  as  that  from  b'  to  o'.  Also,  because  the  tangency  of  the  stable  and  unstable 
manifolds  at  o  is  quadratic,  we  have  that  the  curve  segment  fod  is  shorter  than  the  curve 
segment  coa  by  V  fii-  Because  e  is  small,  this  also  implies  that  the  segment  / o' d'  is 
shorter  than  c'o'a'  by  By  putting  these  facts  together  we  summarize  the  relevant 
information  as  follows: 

(oV)  =  (T?),  (6) 

(c'o'a')  =  (f'o"d"m/4W),  (7) 

where  the  superscripted  bar  denotes  the  length  of  the  segment. 

We  now  wish  to  obtain  an  estimate  of  the  ratio  of  the  measure  of  the 
attractor  contained  in  the  region  d"e"f'o"d"  to  the  measure  of  the  attractor  contained  in 
the  region  a'o'c'b'a'.  To  do  this,  imagine  iterating  the  shaded  region  at  the  tangency  o 
forward  many  iterates  in  time  so  that  o  maps  to  the  point  o +  in  Fig.  5.  The  shaded  region 
has  now  been  greatly  stretched  out  along  the  direction  of  the  unstable  manifold  of  B. 
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Because  dfj.(e)/de  ~  e^1,  we  conclude  that  the  amount  of  attractor  measure  per  unit 

horizontal  length  contained  in  the  shaded  region  emanating  from  o+  is  larger  in 
d'e"f'o"d'  than  in  doc'b'd  by  the  factor  p\~r .  The  same  result  applies  for  other  iterates  of 


the  shaded  region.  Thus  using  Eq.  (7)  we  obtain 


ft 


1/2 


^  R\-y 
Ida'o'c'b'd)  ~  pi  Pi  ’ 


(8) 


where  \i{d'e"f'o"d' )  represents  the  measure  in  d'e"f'o"d'  while  jMa'o'c'b'd)  represents 
the  measure  in  do'c'b'd. 

Because  the  region  d' e"f'o"d‘ ’  maps  to  defod,  the  measure  of  the  attractor 
inside  d'e"f'o"d'  is  just  H(£fh)-  Similarly,  the  measure  of  the  attractor  inside  do'c'b'd  is 
fxie).  Thus,  because  fi{£)~£Y,  we  also  have 


M(dvyy<r)  of 

fj{do'c'b'd)  -Pi- 


(9) 


By  combining  Eqs.  (8)  and  (9)  we  have 

p121=(Pim 

which  yields  the  desired  result,  Eq.  (5)  for  y. 

Finally,  although  we  have  not  done  numerical  experiments  with  noisy 
D -dimensional  maps  (D> 2),  we  conjecture  that  scaling  form  (2)  will  apply  to  many 
situations  in  higher  dimensions  also,  if  Eq.  (1)  applied  in  the  absence  of  noise.  The 
projection  procedure  used  to  reduce  the  two-dimensional  case  to  the  one-dimensional 
case  could  also  be  applied  to  higher  dimensional  systems.  The  main  requirement  for 
scaling  form  (2)  to  apply  to  a  higher  dimensional  situation  is  that  the  measure  of  the 
attractor  accumulate  as  /r(£)~£7  near  the  edge  of  the  attractor.  Reference  5  (1987) 
discussed  the  generalization  of  Eq.  (4)  to  smooth  D-dimensional  maps,  for  the  case  where 
there  is  an  attractor  containing  an  unstable  periodic  orbit  with  D— 1  unstable  eigenvalues 
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and  one  stable  eigenvalue,  and  the  attractor  undergoes  the  D  -dimensional  analog  of  a 
heteroclinic  tangency  crisis.  This  particular  generalization  conformed  to  the  required 
scaling  for  the  measure.  Physical  experiments  have  been  done  producing  chaotic 
attractors  in  the  Poincare  section  having  correlation  dimensions  greater  than  2  which  also 
demonstrated  the  required  scaling  for  the  measure.^  However,  in  those,  and  indeed  most, 
higher  dimensional  cases  there  is  as  yet  no  theoretical  basis  for  predicting  the  critical 
exponent  y.  This  is  one  reason  why  noise  may  be  useful  as  an  experimental  probe  in 
physical  systems,  where  scaling  form  (2)  may  complement  scaling  form  (1)  in 
investigating  critical  phenomena. 
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III.  THE  NOISY  LOGISTIC  MAP 


In  this  section,  we  estimate  analytically  the  characteristic  time  T  of  a 
particular,  one-dimensional  noise-induced  crisis  resulting  in  attractor  destruction.  This 
case  was  first  discussed  in  Ref.  9.  We  show  that  the  dynamics  of  successive  iterates  of  the 
system  near  the  crisis-mediating  fixed  point,  which  contribute  to  high  order  effects  such 
as  orbits  leaving  and  reentering  the  basin  of  attraction,  do  not  affect  the  validity  of  scaling 
form  (2).  We  also  discuss  numerical  experiments  performed  on  the  system,  which 
demonstrate  that:  (i)  the  transient  length  T  is  exponentially  distributed,  as  implied  by  the 
derivation  of  scaling  form  (2),  and  (ii)  our  analytical  estimate  of  f,  which  conforms  to 
scaling  law  (2),  agrees  with  the  numerical  results  over  a  range  of  parameters,  suggesting 
that  the  asymptotic  scaling  form  (2)  is  in  fact  observable  in  practice. 

The  noisy  logistic  map  is  given  by 

xn+i=pxn(l-Xn)  +  Zn.  (10) 

We  consider  the  case  where  the  Zn  are  independent,  identically  distributed  Gaussian 
random  variables  with  variance  <7^.  [In  this  dissertation,  we  will  denote  a  Gaussian 
probability  density  with  mean  m  and  variance  cP-  as  N(»;m,cfi)',  a  uniform  probability 
density  on  an  interval  ( a,b )  will  be  denoted  U(»;a,b).]  It  is  well  known  that  the 
deterministic  (Zn  identically  0)  analog  of  Eq.  (10)  at  has  a  crisis  at  p=Pc=4  resulting  in  the 
destruction  of  a  chaotic  attractor  covering  the  interval  [0,1].  Further,  atp=4  the  chaotic 


attractor  has  invariant  measure  given  by11 


p(x)  = 


1 

Wr(l-r) 


For  p  slightly  less  than  4,  the  attractor  is  reduced  in  size  to  the  interval  [4—p,p/A].  For  our 


purposes,  in  what  follows  we  assume  that  the  measure  for  p  slightly  less  than  4  can  be 

approximated  by 

p(x)  =  j 1  -€Kp/A~x) , 

jt^p/A-x 
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for  x  near  /?/4;  for  <7  small  enough,  we  do  not  care  about  the  form  of  p(x)  outside  this 
neighborhood.  The  dynamical  density  P^+Z(w)  can  be  obtained  by  convolving  P^j)  and 


p(x).  This  density  has  a  nonzero  tail  to  the  right  of  w=l,  i.e.  outside  the  attractor’s  basin 
of  attraction.  We  estimate  a  contribution  R  (escapes  to  the  right)  to  the  overall  escape  rate 
by  integrating  Px+^w) over  obtaining 


(1  Pc-P' 


(ID 


where  Dk(x)  is  the  parabolic  cylinder  function.12  Equation  (11)  was  obtained  in  Ref.  9, 


and  was  used  to  suggest  that  scaling  laws  similar  to  those  applicable  in  noisy  period¬ 
doubling  cascades  and  intermittency^  should  also  apply  to  noise-induced  crises.  Because 
the  critical  exponent  y  for  this  one-dimensional  crisis  is  ^  and  taking  Eq.  (11) 


conforms  to  scaling  law  (2). 

By  applying  the  map  (10)  to  the  dynamical  density,  it  is  possible  to  obtain 
analytical  estimates  for  the  fraction  R  of  the  escapes  contributing  to  R  that  fail  because 
the  noise  reinjects  the  orbit  into  the  basin  of  attraction  on  the  next  iterate,  and  also  for  the 
rate  L  of  escape  from  the  left  side  of  the  attracting  interval.  The  derivation  is  outlined  in 
Fig.  6.  and  done  in  detail  in  Appendix  A.  The  right-hand  escape  rate  R  was  obtained  by 
integrating  the  dynamical  density  Px+Z (w)  over  [1»°°)  Fi§-  6(a)].  By  applying  the 
deterministic  part  of  the  map  (10)  to  the  density  Px+Z(w)’  one  °htains  the  density  of 
orbits  before  they  are  kicked  by  noise  on  the  next  iterate,  P x+z(w)  tc^*  Fig.  6(b)]. 
Convolving  the  noise  density  with  that  part  of  Px+Z(w)  to  fifiht  °f  ancl 
integrating  the  convolution  over  (— < °°,0],  one  obtains  an  estimate  of  the  left-hand  escape 
rate  L  [cf.  Fig.  6(c)].  This  approach  assumes  that  most  of  the  left-hand  escapes  will  in  fact 
be  not-quite-successful  right-hand  escapes.  Because  of  the  expansion  produced  by  the 
unstable  fixed  point  at  0,  it  becomes  progressively  less  likely  for  an  orbit  to  be  kicked  out 
of  the  left  side  of  the  attractor  on  successive  iterates.  Thus  the  assumption  seems 
warranted.  Alternatively,  convolving  the  noise  density  with  that  part  of  P x+z(w) t0  the 
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left  of  w=0,  and  integrating  the  convolution  over  (0,°o]  yields  an  estimate  of  the  rate  R  at 
which  right-hand  escaping  orbits  are  kicked  by  the  noise  back  into  the  basin  of  attraction 
(cf.  Fig.  6(d)). 


(a) 


P x+z(w) 


(b) 


FIG.  6.  (a)  The  right-hand  escape  rate  for  the  noisy  logistic  map  is  obtained  by 
integrating  the  dynamical  density  />y+z(w)  over  [1,°°).  (b)  Applying  the  deterministic 

part  of  map  (1)  yields  the  density  F^+Z(w)  of  orbits  before  the  next  noise  kick,  (c)  and 
(d)  The  mapped  density  selectively  convolved  with  the  noise  density  allows 

estimation  of  the  left-hand  escape  rate  L  and  the  right-hand  failed-escape  rate  R . 
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These  higher-order  contributions  to  the  overall  escape  rate  can  be  expressed  as 
asymptotic  expansions  (see  Appendix  A): 


0-1/2 
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Therefore,  an  estimate  of  the  characteristic  time  t  can  be  written  as 

r=[R-~R  +L]~l,  (14) 

where  one  can  see  from  Eqs.  (11)— (13)  that  scaling  form  (2)  with  7=  2  aPP^es-  l11 

addition  to  the  three  terms  R,  R ,  and  L,  there  will  also  be  corrections  due  to  higher  order 
effects  such  as  multiple  crossings  of  the  fixed  point  on  successive  iterations  of  the  noisy 
map. 

The  derivation  of  the  correction  terms  R  and  L  shows  that  for  a  noise 
density  Pz(z)  that  vanishes  fast  enough  with  z,  the  dynamics  around  the  fixed  point, 

which  only  linearly  distort  the  overall  picture  from  iterate  to  iterate,  do  not  affect  the  final 
scaling  form  because  we  only  assume  a  linear  relation  between  Pc~P  and  the  distance 
from  the  attractor  to  the  mediating  fixed  point;  different  constants  of  proportionality  at 
each  iterate  generate  contributions  to  t  that  all  scale  like  Eq.  (2).  Thus  the  higher-order 
terms  do  not  change  the  scaling  form  (2). 

Note  also  that  although  R  is  usually  small  compared  to  the  other  terms  in 
Eq.  (14),  L  tends  to  dominate  R  as  (pc  -  p)/o  increases  [R~L  when  (pc-p)/o=  12].  This 
result  differs  from  a  conclusion  of  Ref.  9,  where  an  estimate  of  L  based  on  numerical 
experiments  was  given  that  underestimated  its  relative  importance  for  large  (pc  —p)/o. 
Thus,  the  noise  “causes”  the  crisis  in  different  parts  of  the  phase  space,  depending  on  its 
strength.  This  property  could  be  important  in  the  investigation  of  higher  dimensional 
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physical  systems,  where  the  boundaries  of  basins  of  attraction  can  be  much  more 
complicated. 

Numerical  experiments  were  also  performed  with  the  noisy  logistic  map  in 
the  following  way.  Initial  conditions  were  chosen  at  random  in  the  deterministic 
attractor’s  basin  of  attraction,  and  were  preiterated  using  the  deterministic  map  at  a 
specified  parameter  value,  to  obtain  the  natural  measure  of  the  attractor.  Then  the  noisy 


map  was  applied,  using  a  specified  noise  variance,  until  the  orbit  escaped  from  the  basin 
of  attraction.  This  was  determined  by  waiting  for  the  orbit  to  reach  the  region  x<-l, 
where,  for  the  values  of  <7  used,  the  chances  of  the  noise  reinjecting  the  orbit  into  the 
interval  [0,1]  were  negligible.  Various  diagnostics  allowed  determination  of  the  relative 
contributions  of  the  various  terms  in  Eq.  (14).  An  ensemble  of  transients  was  sampled  for 
each  set  of  parameters,  allowing  statistical  estimation  of  the  characteristic  time  f.  This 
general  procedure  was  used  for  all  of  the  numerical  results  presented  in  this  dissertation, 
with  specific  criteria  for  deciding  the  end  of  a  transient  response  depending  on  the 
particular  system.  Use  of  several  different  end  criteria  in  a  given  system  only  changed  the 
constant  of  proportionality  in  scaling  laws.  For  all  of  the  numerical  experiments, 
pseudorandom  numbers  satisfied  distribution  and  independence  requirements  to  high 
statistical  confidence  levels. 

Figure  7  shows  a  typical  empirical  distribution  of  transient  lengths  for  the 

t 

noisy  logistic  map.  The  cumulative  distribution  function  Fj(t)=  \PT(s)ds  for  transient 

0 


lengths  T  theoretically  should  be  of  the  form 
Fj{t)  =  1  -  exp(-t/x),  t> 0, 

where  T  is  the  characteristic  time  at  the  specified  parameter  value  and  noise  strength.  The 
sampled  distribution  is  compared  with  this  functional  form,  using  the  sample  mean  <t>  as 
an  unbiased  estimator  of  t. 
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FIG.  7:  The  time  (in  iterations)  taken  by  trajectories  of  the  noisy  logistic  map  to 
escape  (specifically,  to  reach  the  region  x<-l)  from  the  deterministic  basin  of 
attraction  is  exponentially  distributed.  The  straight  line  shows  the  theoretical  form  of 
the  distribution,  which  is  statistically  consistent  with  the  data.  A  Kolmogorov- 
Smimov  test  indicated  that  the  probability  of  the  deviations  from  the  theoretical 
distribution  being  as  large  as  indicated  due  to  chance  is  0.38.  In  this  and  later  figures, 
error  bars  indicate  an  approximately  68%  probable  (one  standard  deviation  of  the 
mean)  confidence  interval  for  the  estimated  quantity. 
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In  Fig.  7  and  what  follows,  we  denote  estimators  of  statistical  quantities  by  a  circumflex. 
The  comparison  is  statistically  significant,  and  combined  with  the  results  from  63  other 
ensembles  at  different  parameter  values,  support  the  conclusion  that  the  characteristic 
transient  lengths  T  are  exponentially  distributed  for  noise-induced  crises.  This  tends  to 
vindicate  the  independence  assumptions  used  in  obtaining  scaling  law  (2). 

The  numerical  results  are  summarized  in  Fig.  8.  Rewriting  scaling  law  (2), 


we  have 


%af~z 


(20 


which  implies  that  as  (pc-p)/c  is  varied,  to?  should  fall  on  a  single  curve  [the  graph  of 
g(*)]  regardless  of  the  value  of  o.  The  results  in  Fig.  8  conform  to  this  prediction  for  a 
range  of  c  spanning  several  orders  of  magnitude.  Moreover,  the  numerical  results  agree 
with  the  analytically  predicted  form  of  #(•)  based  on  Eqs.  (11) — (14).  Although  not  shown, 
the  individual  contributions  R,  R  ,  and  L  separately  agree  with  the  corresponding 
numerical  results.  For  the  results  in  Fig.  8,  kmax-3  was  used  to  evaluate  Eqs.  (11)  and 
(12);  R  was  about  0.1/?,  and  L  ranged  from  about  0.1/?  at  (pc  —  p)/ <7=0  to  about  1.2/?  at 


(pc-p)i<y=  14- 
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FIG.  8:  The  results  of  the  numerical  experiments  with  the  noisy  logistic  map  are 
statistically  consistent  with  the  theoretically-predicted  scaling  function  given  in  Eqs. 
(11)— (14),  for  a  three  order-of-magnitude  variation  in  the  strength  of  the  additive 
noise. 
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IV.  NUMERICAL  EXAMPLES  USING  NOISY  TWO-DIMENSIONAL  MAPS 


In  this  section,  we  illustrate  with  numerical  examples  the  application  of 
scaling  law  (2)  to  two  two-dimensional  dissipative  maps  with  noise,  and  demonstrate  a 
restricted  form  of  (2)  reminiscent  of  scaling  law  (1),  where  the  noise  strength  takes  the 
role  of  the  parameter  excess  p-pc> 0. 


A.  Henon  map;  heteroclinic  tangency  crisis 

The  Henon  map  was  modified  with  additive  noise  to  produce  the  system 
2  r  m 
Xn+i  —  p  —  Xjj  —  Jyn  +  Zyi 


y/l+1  Xfl 

where  the  Z n=(Z^\  Z  ®)  are  independent,  identically  distributed  two-vectors. 


For  J=  -0.3,  p=pc~  1.426  921  114,  the  deterministic  version  of  Eq.  (15)  has  a  heteroclinic 
tangency  crisis,  caused  by  an  unstable  fixed  point  [corresponding  to  B  in  Fig.  3(a)],  which 
results  in  destruction  of  the  chaotic  attractor  existing  for  p<pc.  Based  on  the  eigenvalues 
of  an  unstable  fixed  point  on  the  attractor  (corresponding  to  A  in  Fig.  3(a)),  Eq.  (4)  yields 
a  critical  exponent  for  the  deterministic  scaling  law  (1)  of  ><=0.8557.  Previous  numerical 
experiments  with  the  deterministic  system  agreed  with  these  theoretical  results.5  This 
value  of  y  was  used  to  scale  the  data  from  numerical  experiments  on  the  noisy  system  of 

Eq.  (15)  for  comparison  with  scaling  law  (2'),  with  the  results  shown  in  Fig.  9.  For  these 
experiments,  the  noise  density  used  was  Pj(x,y)=N(x;0,cP-)S(y).  We  are  no  longer  able  to 

provide  a  theoretical  determination  of  the  function  g(*)  (as  was  done  in  Chapter  III)  for 
comparison  with  the  data,  but  the  statistical  consistency  of  the  data  over  three  orders  of 
magnitude  for  the  noise  strength  tends  to  confirm  the  applicability  of  Eq.  (2)  to  this 
system. 
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FIG.  9.  Scaling  form  (2')  is  observed  for  a  three  order-of-magnitude  range  of  noise 
strength  in  the  noisy  Henon  map  operated  near  a  heteroclinic  tangency  crisis. 
y~  0.8557  based  on  the  eigenvalues  of  the  period-1  saddle  on  the  attractor  and  Eq. 
(4). 
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The  non-universality  of  the  function  g(*)  in  Eq.  (2)  is  demonstrated  in  Fig. 
10,  where  the  results  of  numerical  experiments  using  several  different  noise  density 
functions  (see  figure  caption)  are  compared.  Although  varying  <7  leaves  data  from 
experiments  with  the  same  noise  density  along  a  single  curve,  the  curves  produced  by  the 
different  densities  have  no  apparent  relation  to  one  another.  One  way  of  addressing  the 
lack  of  a  theoretical  prediction  for  #(•)  in  experiments  is  to  operate  the  system  at  the 
parameter  value  corresponding  to  the  deterministic  crisis.  Then  Eq.  (2)  is  restricted  to  a 
form  akin  to  the  deterministic  scaling  law  (1),  with  the  noise  strength  staking  the  role  of 
the  parameter  excess  p-pc' 

T~<j-rg(0)~<j-y.  (2") 

A  comparison  of  the  numerical  results  from  the  noisy  Henon  map  with  the  restricted 
scaling  (2")  is  given  in  Fig.  11,  where  the  statistical  consistency  again  supports  the 
applicability  of  Eq.  (2)  to  this  system. 


29 


0.0  1.0  2.0  3.0  4.0 

Pc~P 

a 


FIG.  10:  Different  noise  densities  result  in  different  functions  g(*)  in  Eq.  (2),  even  for 
the  same  system,  in  this  case  the  noisy  Henon  map  near  the  heteroclinic  tangency 
crisis  dealt  with  in  Fig.  9.  The  three  types  of  symbol  each  relate  to  a  different  form  of 
the  density  of  the  additive  noise;  the  open  and  closed  symbols  of  a  given  type 
demonstrate  that  the  scaling  (2)  pertains  over  a  range  of  c  within  each  type  of  noise. 
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FIG.  1 1.  Restricting  p=pc  fixes  the  argument  of  g(a)  in  Eq.  (2),  resulting  in  a  scaling 
similar  to  Eq.  (1).  The  slope  of  the  straight  line  is  -7= -0.8557  based  on  the 
eigenvalues  of  the  period- 1  saddle  on  the  attractor  and  Eq.  (4). 


B.  Ikeda-Hammel-Jones-Maloney  map;  homoclinic  tangency  crisis 

The  Ikeda-Hammel-Jones-Maloney  map,  a  model  for  a  nonlinear  optical 
ring  resonator13  was  modified  by  additive  noise  to  produce  the  system 

wn+ 1  =  A  +  Bwnexp[iK-ia/(l+\wn\2)]  +  Zn,  (16) 

where  wn=xn+iyn;  x,  y.  A,  B,  K,  a,  Zn  are  all  real;  and  Zn  are  independent,  identically 
distributed  random  variables  with  a  uniform  density  U(x;-<J,cf).  For  A=AC*=1.0024, 5=0.9, 
>o=0.4,  and  or=6.0,  the  deterministic  system  corresponding  to  map  (16)  has  a  homoclinic 
tangency  crisis  mediated  by  an  unstable  fixed  point  (corresponding  to  B  in  Fig.  3(b)), 
which  results  in  destruction  of  the  attractor  existing  for  A<AC.  (Note  that  in  this  case,  A 
takes  the  role  of  the  generic  crisis  parameter  p.)  Based  on  the  eigenvalues  of  the  unstable 
fixed  point,  Eq.  (5)  yields  a  critical  exponent  for  the  deterministic  scaling  law  (1)  of 
yts 1.705.  The  relatively  large  value  of  the  exponent  y  limited  the  dynamic  range  of  noise 
that  could  be  investigated  numerically,  due  to  the  extremely  long  transients  obtained  for 
smaller  values  of  <r.  (In  fact,  to  obtain  even  somewhat  less  than  three  orders  of  magnitude 
variation  in  <7 ,  censored  sampling  techniques  were  necessary.  In  censored  sampling  to 
determine  a  characteristic  time,  an  ensemble  of  m  orbits  is  advanced  in  parallel  iterate  by 
iterate.  The  mean  lifetime  is  deduced  statistically14  after  only  the  first  r<m  orbits  end 
their  transient  responses,  resulting  in  a  large  computational  saving  relative  to  allowing  all 
of  the  ensemble  of  r  orbits  to  escape.  The  resulting  statistical  confidence  is  equivalent  to 
the  r  out  of  r  case.  In  this  and  the  following  examples,  typically  r=128,  m= 256.) 
Nevertheless,  the  results  of  the  numerical  experiments,  when  compared  to  restricted 
scaling  form  (2")  as  in  Fig.  12,  suggest  that  this  system  also  exhibits  the  predicted 
scaling. 
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FIG.  12.  The  restricted  scaling  law  (2")  is  consistent  with  the  results  of  this  numerical 
experiment  on  the  noisy  Ikeda-Hammel-Jones-Maloney  map  near  a  homoclinic 
tangency  crisis.  The  slope  of  the  straight  line  is  —  y~  —1.705  based  on  the  eigenvalues 
of  the  period- 1  saddle  mediating  the  crisis  and  Eq.(5). 


V.  POTENTIAL  COMPLICATIONS  TO  OBSERVING  SCALING  LAW  (2) 

Although  scaling  law  (2)  agreed  with  the  results  of  the  numerical 
experiments  on  the  systems  discussed  in  the  previous  chapter,  we  have  also  found 
examples  where  agreement  with  the  scaling  is  less  immediately  apparent.  In  this  section, 
we  discuss  two  such  examples  that  illustrate  problems  that  might  be  encountered  in 
physical  experiments  with  noisy  systems.  The  first  example  (Chapter  V  A)  shows  how 
frequently-encountered  properties  of  the  measure  of  chaotic  attractors,  that  were  not 
explicitly  addressed  in  the  derivation  in  Chapter  II,  can  lead  to  results  that  cluster  around, 
rather  than  on,  the  curve  g(»)  in  Eq.  (2).  The  second  example  (Chapter  V  B)  shows  that 
the  range  of  noise  strengths  over  which  scaling  law  (2)  can  be  observed  may  be  restricted. 
This  happens  when  there  are  mechanisms,  different  from  the  tangencies  of  the 
deterministic  crisis,  by  which  the  noise  can  cause  orbits  to  escape. 

A.  Ragged  measure  of  deterministic  attractor 

The  analysis  of  Chapter  II B  leading  to  Eq.  (5)  for  a  homoclinic  tangency 
crisis  showed  that  the  striations  of  the  strange  attractor  (which  has  a  Cantor-set  structure) 
near  its  tangencies  with  the  stable  manifold  of  the  periodic  point  mediating  the  crisis, 
accumulate  at  the  geometric  rate  /%,  the  stable  eigenvalue  of  the  periodic  point  B  of  Fig. 
3(b).  This  means  that  as  one  considers  the  attractor  measure  as  a  function  of  distance 
away  from  one  of  the  tangencies,  measured  perpendicularly  relative  to  the  stable 
manifold,  it  has  a  “stair-step”  quality  with  successive  steps  being  l/p2  wider  than  the  last. 
This  effect  has  been  observed  in  numerical  experiments  on  deterministic  crises.  In 
particular,  in  log-log  plots  of  T  versus  (p~Pc )>  oscillations  periodic  in  log\o(p-pc)  were 
seen  to  be  superimposed  on  the  logarithmically  linear  behavior  of  scaling  law  (1),  where 
the  period  of  the  oscillation  is5  logioiHfii).  For  cases  where  /J2«l»  this  effect  is 
particularly  pronounced,  and  we  refer  to  it  as  resulting  from  a  ragged  measure.  In  noisy 
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systems,  a  ragged  measure  can  also  complicate  the  observation  of  scaling  law  (2)  as  we 
now  demonstrate  with  another  numerical  example  based  on  the  noisy  Henon  map,  Eq. 
(15). 

For  7=0.3,  p=pc=2.124  672  450  the  deterministic  analog  of  Eq.  (15) 
undergoes  a  homoclinic  tangency  crisis  mediated  by  an  unstable  period-3  orbit,  which 
results  in  destruction  of  the  chaotic  attractor  existing  for  p^Pc-  The  eigenvalues  of  the 
mediating  saddle  are  j8i=10.9  and  £2=0.00247;  the  small  stable  eigenvalue  leads  to  a 
pronounced  periodic  component  in  numerical  results  plotted  as  log io('T)  versus 
log \oip~~Pc\  with  period  fogio(l//fc)~2-6  (base  ten)  log  units. ^  An  analogous  periodicity 
is  noted  in  the  results  of  numerical  experiments  intended  to  test  the  restricted  scaling 
form  (2")  on  the  noisy  system,  as  shown  in  Fig.  13  [where  the  noise  was  U(x;—o,cf)].  The 
oscillations  superimposed  on  the  straight-line  scaling  implied  by  Eq.  ( 2 ")  have  an 
amplitude  we  denote  Alog\o(t). 

The  observation  of  oscillations  in  the  the  noisy  system  probably  depends 
on  more  than  the  presence  of  a  ragged  measure.  In  the  deterministic  case,  a  parameter 
excess,  p-pc> 0,  leads  to  fingers  of  the  attractor  remnant  poking  over  the  stable  manifold 
(like  region  ab  in  Fig.  4  for  the  heteroclinic  case)  where  each  finger  is  the  image  under 
the  map  of  the  previous  finger.  In  the  deterministic  case,  it  suffices  to  look  at  the  measure 
in  a  single  finger,  because  once  an  orbit  crosses  the  stable  manifold,  it  stays  across  the 
stable  manifold  on  successive  iterates.  Therefore,  as  a  region  in  one  finger  with  a  large 
value  of  dfi(e)/d£  [corresponding  to  a  region  of  higher-than-average  slope  on  a  plot  of 
log  io(f)  versus  logio(p-pc)]  just  crosses  the  stable  manifold,  regions  of  correspondingly 
high  dp(e)/de  in  all  the  other  fingers  will  simultaneously  cross  the  stable  manifold  as 
well.  However,  the  situation  is  different  in  the  noisy  case.  For  a  given  noise  strength  a 
orbits  can  escape  if  they  are  located  within  a  distance  of  order  cr  from  the  stable  manifold 
of  periodic  orbit  B. 
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FIG.  13.  Results  of  the  numerical  experiment  on  the  noisy  Henon  map  near  a 
homoclinic  tangency  crisis  mediated  by  a  period-3  saddle.  The  restricted  scaling  (2  0 
shows  evidence  of  the  raggedness  of  the  attractor’s  natural  measure.  logio(t)  has  a 
component  periodic  in  log\o(l/fi2)  =  /ogio(404.9...)  =  2.61  with  an  amplitude  of 
Alogioix).  The  slope  of  the  straight  lines  are  -y  ~  -0.83  based  on  the  eigenvalues  of 
the  period-3  saddle.  The  plot  symbols  at  upper  right  indicate  the  values  of  <7 for  which 
the  corresponding  curves  of  Fig.  14  were  computed. 
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Each  finger  (now  analogous  to  region  aoc  in  Fig.  5)  represents  an  independent 
opportunity  for  escape.  Furthermore,  the  region  in  each  finger  which  is  accessible  for 
escape  is  not  simply  the  image  of  the  corresponding  region  in  the  previous  finger.  Each 
finger  makes  a  contribution  to  the  overall  escape  rate  that  is  periodic  in  logio(o),  but  the 
phases  of  these  separate  oscillations  are  different.  This  would  tend  to  reduce  the 
amplitude  of  any  observable  oscillations.  However,  if  the  measure  in  one  finger 
dominates  that  in  the  others,  oscillations  such  as  those  in  Fig.  13  can  be  seen  regardless  of 
the  phase  relationships.  Such  a  situation  seems  relatively  common. 

Oscillations  superimposed  on  the  restricted  scaling  (2  )  produce  even 
stronger  effects  on  results  relevant  to  the  general  scaling  (2).  Consider  the  effects  of  the 
ragged  measure  in  this  example  for  o  held  constant  as  pc  —  p  is  increased  from  zero.  The 
quantity  io?  increases  as  g(»)  increases.  The  starting  point  of  the  curve,  the  value  of  g( 0), 
should  be  obtainable  from  the  data  in  Fig.  13.  However,  note  that  different  values  of  g(0) 
are  obtained  for  different  <r  due  to  the  oscillations  about  scaling  law  (2  ).  Thus,  for 
different  values  of  <J  the  quantity  t gY  should  follow  different,  parallel  curves  as  (pc  ~P)/& 
is  varied.  This  is  essentially  what  is  seen  in  Fig.  14,  where  the  results  fall  in  a  band  of 
approximate  width  Alog  the  same  as  the  amplitude  of  the  superimposed  oscillations 
of  Fig.  13.  In  general,  results  for  <r  corresponding  to  a  trough  in  Fig.  13  fall  in  the  lower 
part  of  the  band  in  Fig.  14,  while  results  for  a  corresponding  to  a  crest  in  Fig.  13  fall  in 
the  upper  part  of  the  band  in  Fig.  14.  Obviously,  if  extensive  data,  as  in  Fig.  14,  are  not 
taken  from  a  system  with  the  properties  of  this  example,  one  would  have  difficulty 
confirming  the  applicability  of  Eq.  (2).  This  could  also  result  in  an  erroneous 
determination  of  the  critical  exponent  y. 
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FIG.  14.  The  scaling  (2)  is  obeyed  only  approximately  for  the  noisy  Henon  map  near 
the  homoclinic  tangency  crisis  dealt  with  in  Fig.  13.  The  results  fall  in  a  band  of 
width  Alog\o(x),  (the  same  as  the  amplitude  of  the  wiggles  in  Fig.  13),  rather  than  on 
a  single  curve. 


B.  Competing  routes  of  escape 

Another  potentially  common  obstacle  to  observing  scaling  behavior  (2)  in 
physical  systems  is  the  presence  of  alternative  routes  for  orbits  to  escape  under  the 
influence  of  noise,  i.e.,  routes  other  than  the  crisis-producing  tangencies.  We 
demonstrate  this  with  another  crisis  of  the  Ikeda-Hammel-Jones-Maloney  map  given  by 
Eq.  (16).  For  A=0.85,  B= 0.9,  k=0A,  and  a=ac~ 7.268  849,  the  deterministic  system 
corresponding  to  map  (16)  has  a  heteroclinic  tangency  crisis  mediated  by  an  unstable 
period-5  orbit,  which  results  in  crisis-induced  intermittency  known^  to  conform  to  Eq.  (1) 
for  a>ac  (in  this  case  a,  rather  than  A,  takes  the  place  of  the  generic  crisis  parameter  p). 
The  results  of  numerical  experiments  intended  to  confirm  the  applicability  of  restricted 
scaling  law  ( 2 ")  for  the  noisy  map  [the  noise  was  again  U(x;-a,&)\  are  shown  in  Fig.  15. 
The  fit  to  the  predicted  exponent  y  is  only  good  for  very  small  <7  (<7<10“3). 

The  reason  for  this  failure  is  shown  schematically  in  Fig.  16.  The  stable 
manifold  colliding  with  the  chaotic  attractor  makes  a  number  of  close  approaches  to  the 
chaotic  attractor  that  nevertheless  do  not  become  tangencies  for  p=Pc •  These  close 
approaches,  although  having  no  effect  on  the  observability  of  scaling  law  (1)  for  the 
deterministic  system,  provide  an  additional  route  of  escape  for  the  noisy  trajectories,  at 
larger  values  of  <7,  which  was  not  accounted  for  in  deriving  Eq.  (2).  Therefore,  at  large 
values  of  <7  the  observed  rate  of  escape  has  two  contributions,  one  from  the  crisis- 
producing  tangencies  and  one  from  the  close  approaches  of  the  stable  manifold,  resulting 
in  a  smaller-than-expected  value  of  t.  For  small  values  of  <7,  the  close  approaches  do  not 
contribute  to  the  rate  of  escape,  and  f  conforms  to  Eq.  (2").  Of  course,  it  is  not  at  all 
surprising  that  a  given  system  will  have  an  upper  bound  on  <7  for  Eq.  (2)  to  apply,  but 
competing  routes  of  escape  make  that  bound  more  stringent  than  it  would  otherwise  be, 
and  provide  a  potential  pitfall  in  the  experimental  determination  of  critical  exponents  of 
noise-induced  crises. 
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FIG.  15.  The  restricted  scaling  (2")  is  violated  for  large  a  ( i.e 
a>~10~3)  in  this  numerical  experiment  with  the  noisy  Ikeda-Hammel-Jones-Maloney 
map  near  a  heteroclinic  tangency  crisis  mediated  by  a  period-5  saddle.  The  slope  of 
the  straight  line  is  -y~  -1.236  based  on  the  eigenvalues  of  the  period-5  saddle  on  the 
attractor  and  Eq.  (4). 
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FIG.  16.  Schematic  of  a  “near  tangency”  unrelated  to  the  crisis-producing  tangency. 
The  arrows  indicate  a  route  for  noisy  trajectories  to  escape  that  is  not  accounted  for  in 
the  derivation  of  Eq.  (2).  Thus,  for  large  a  an  additional  rate  of  escape  will  decrease  f 
from  what  is  predicted  on  the  basis  of  the  crisis  alone. 
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VI.  APPLICABILITY  OF  SCALING  LAW  (2)  TO  FLOWS 

Thus  far,  we  have  been  rather  arbitrary  in  adding  noise  to  the  example 
systems,  and  have  ignored  how  noise  terms  might  arise  physically.  In  particular,  the  noise 
densities  have  been  assumed  to  be  independent  of  phase-space  position.  Although  it  is 
possible  to  imagine  physical  systems  for  which  a  discrete  map  is  explicitly  derivable  from 
first  principles  (the  Ikeda-Hammel-Jones-Maloney  map  is  one  example),  and  where  noise 
might  appear  in  this  relatively  benign  fashion,  in  general  we  expect  a  more  complicated 
situation.  Indeed,  much  of  the  justification  for  studying  discrete  maps  is  that  they  can  be 
viewed  as  occurring  in  a  Poincare  section  of  the  higher-dimensional  phase  space  of  a 
flow.  However,  noisy  flows  modeling  physical  systems  will  in  general  not  reduce,  in 
Poincare  sections,  to  noisy  maps  of  the  type  seen  in  the  examples  so  far.  In  particular,  the 
noise  applied  to  such  a  model  to  model  the  physical  system  will,  in  general,  not  be  simply 
additive  nor  will  it  be  independent  of  the  coordinates  in  the  surface  of  section.  In  this 
chapter  we  present  an  example  of  a  noisy  flow  that  is  physically  well-motivated,  and  use 
it  to  suggest  that  Eq.  (2)  is  also  applicable  to  the  more  realistic  situation  where  the  density 
function  of  the  noise  in  a  noisy  map  depends  on  phase- space  position. 

We  consider  a  noisy  version  of  the  Duffing  oscillator,  a  model  for  the  one¬ 
dimensional  motion  of  a  point  particle  in  a  potential  V(x),  which  is  subjected  to  friction, 
an  external  sinusoidal  force,  and  random  forcing  like  that  producing  Brownian  motion. 
This  situation  is  described  by  the  stochastic  differential  equation 

c fix/dt 2  =  -vdx/dt-dV/dx  +  p  sin(cot )  +  £(r),  (I7) 

where  £(f)  is  Gaussian  white  noise  with  intensity  <7,  i.e.,  (£(r)£( 0)  =  cPS(t-t').  We  take 
V(x)  =  aX*/ 4-  -  fix2/!,  v=l ,  a=100,  >3=10,  <y=3.5  and  p  near  pc=0.845.  The  potential  well 
has  two  minima  at  x  =  iV/^a.  For  p<pc,  the  deterministic  analog  of  Eq.  (17)  has  two 
symmetrically  disposed  chaotic  attractors,  one  confined  to  the  well  in  x<0  and  the  other  to 
the  well  in  x>0.  For  p>pc,  there  is  one  chaotic  attractor,  and  orbits  intermittently  switch 
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from  one  well  to  the  other.  ^  When  considered  in  the  t=2n7t/co  Poincare  section,  the 
deterministic  system  undergoes  a  homoclinic  tangency  crisis  mediated  by  an  unstable 
period-3  orbit,  which  results  in  attractor  merging.  Numerical  experiments  confirm  that  the 
crisis  scaling  (1)  applies  to  this  system.5^9**6) 

One  of  the  deterministic  attractors  for  p<pc  is  shown  in  Fig  17(a).  When 
the  stochastic  system  of  Eq.  (17)  is  integrated  as  discussed  in  Appendix  B,  the  envelope 
of  the  orbit  spreads  out  under  the  influence  of  the  noise,  as  shown  in  Fig.  17(b).  Note  that 
because  the  noise  does  not  appear  in  the  argument  of  the  sine  in  Eq.  (17),  the 
stroboscopic  section  is  still  a  valid  Poincare  section  for  the  stochastic  flow.  The 
difference  between  this  example  and  previous  ones  in  this  dissertation  (which  also 
produced  “fuzzy”  attractors)  is  best  shown  by  considering  the  conditional  density  of 
trajectories  under  one  period  of  the  forcing,  given  a  specified  initial  condition,  for  many 
different  realizations  of  the  noise;  two  such  conditional  densities  are  illustrated  in  Figs. 
17  (c)  and  (d).  We  see  that  the  effect  of  the  noise  §(r)  is  different  in  different  parts  of  the 
phase  space,  depending  on  the  local  stretching  and  dissipation  given  by  the  deterministic 
terms  in  Eq.  (17),  resulting  in  a  spatially-dependent  density  for  the  noise  in  the  Poincare 
section.  The  conditional  density  in  the  Poincare  section  tends  to  extend  along  the  unstable 
manifold  of  the  deterministic  attractor.  Although  the  exact  shape  of  the  density  is  phase- 
space-dependent,  the  characteristic  spread  of  the  noise  in  the  Poincare  section  is  of  the 
same  order  over  the  region  of  the  deterministic  attractor,  and  that  characteristic  spread 
scales  linearly  with  the  noise  intensity  parameter  <r.  Therefore,  it  seems  reasonable  to  try 
to  apply  scaling  law  (2)  to  this  and  other  continuous-time  systems. 
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J7S*  l?- (a)  of  the  two  p<pc  attractors  of  the  deterministic  differential  equation 
17),  shown  in  the  t=2nnlco  Poincare  section.  5000  iterates  of  a  trajectory  areXS 
(b)  The  noisy  attractor  for  <7=10-5.  Again,  5000  iterates  were  computed,  but  a  number 
totheXh8 (betWCen  the X<°  and X>0 P°tential  wells  have  already  taken  place,  due 
de*sitieXfttn^anJeter  SfCe)  h0?0Clmic  tangency  crisis-  (c)  (d)  Conditional 

stochastic iffetS 17  ^  °nC  fOTCing  Peri°d’  given  an  initial  condition.  The 
stochastic  differential  equation  was  integrated  for  one  forcing  period  for  800  different 

oise  rea  izations,  starting  from  the  initial  condition  shown  by  the  circle  The 

e=^rCt0ry,1S  Shfwn  by  the  Cross-  ^  noise  density  in  the  Poincare 
section  is  clearly  space-dependent. 
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Figures  18  and  19  show  that  scaling  law  (2),  with  y  determined  from  the 
eigenvalues  of  the  period- 3  orbit,  does  indeed  apply  to  the  stochastic  differential  equation 
(17),  albeit  with  the  complications  discussed  in  Chapter  V  A,  which  are  unrelated  to  the 
fact  that  the  system  is  a  flow,  rather  than  a  map.  Here,  the  scale  of  the  oscillations 
superimposed  on  restricted  scaling  form  (2")  is  particularly  large,  because  the  stable 
eigenvalue  of  the  mediating  period-3  saddle  in  the  Poincare  section  is  so  small 
(02  =  0.000512). 
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FIG.  18.  The  restricted  scaling  (2")  applies  to  the  stochastic  differential  equation  (17), 
with  the  same  type  of  measure-dependent  structure  seen  in  Fig.  13.  The  scale  of 
variation  is  larger  here  because  the  stable  eigenvalue  of  the  mediating  period-3  saddle 
in  the  Poincare  section  is  so  small  (/%  =  0.000512).  The  slope  of  the  straight  line  is 
— y~  -0.703  based  on  the  eigenvalues  of  the  period-3  saddle. 


FIG.  19.  The  scaling  (2')  applied  to  the  results  of  the  numerical  experiment  with  the 
stochastic  differential  equation  (17).  The  scaling  curve  g(*)  is  broadened  as  in  Fig.  14. 


VH.  IMPLICATIONS  FOR  EXPERIMENTS 


In  this  chapter,  we  summarize  the  implications  of  scaling  law  (2)  for 
experiments  on  physical  systems  accompanied  by  noise.  These  implications  are  intrinsic 
to  noisy  systems,  and  are  not  merely  special  (although  common)  problems  like  those 
discussed  in  Chapter  V.  Figure  20  shows  in  three-dimensional  form  the  qualitative 
dependence  of  the  characteristic  time  x  of  noise-induced  crises  on  the  strength  of  the 
noise  <rand  the  parameter  p.  Note  that  the  region  p>pc  is  also  included  in  the  figure.  The 
derivation  in  Chapter  II  did  not  explicitly  address  the  case  of  p>pc,  but  our  experience 
with  numerical  experiments  confirms  the  reasonable  idea  that  for  noisy  systems  with 
( Pc~p)l<y  small  and  negative,  the  behavior  is  continuous  with  that  for 
ipc  -  p)/ g  small  and  positive;  for  ( pc  -  p)/o  large  and  negative,  the  system  behaves  as  in  a 
deterministic  crisis. 

A  typical  experiment  is  to  determine  the  parameter  value  at  which  a  crisis 
occurs,  by  measuring  the  length  of  the  transient  as  a  function  of  the  parameter.  Ideally, 
one  should  do  this  in  a  vertical  plane  <7=0  of  Fig.  20,  where  x  would  go  asymptotically  to 
infinity  at  p=pc.  However,  the  presence  of  noise  limits  the  experiment,  at  best,  to  some 
vertical  plane  <7/0,  where  large  values  of  x  are  found  for  p<pc.  Thus,  there  is  a  tendency 
for  pc  to  be  underestimated  in  systems  with  noise,  and  for  the  system  to  appear  to  disobey 
the  deterministic  scaling  law  (1)  for  p-pc  varied  over  a  wide  dynamic  range. 

If  an  underestimate  of  pc  from  experiments  on  a  supposedly  deterministic 
system  is  used  in  a  subsequent  experiment  where  additional  noise  is  intentionally  added, 
as  a  probe,  to  test  the  applicability  of  scaling  law  (2),  spurious  results  will  again  be 
possible.  For  example,  a  test  of  Eq.  (2")  should  ideally  take  place  in  the  vertical  plane 
Pc  ~  P  -  0  of  Fig.  20,  where  intersection  with  the  surface  yields  a  straight  line  for  ln(x) 
versus  ln(c).  The  previous  underestimate  of  pc  will  lead  to  the  experiment  being  done  in 
a  vertical  plane pc-p>0,  where  the  ln(x)  versus  relation  ln(<7)  will  be  nonlinear. 
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FIG.  20.  The  surface  displays  the  parameter  and  noise  dependence  of  the  transient 
length  for  a  notional  scaling  function  g(»).  The  thick  line  gives  the  restricted  scaling 
(2").  Curves  in  vertical  planes  parallel  to  the  thick  line  correspond  to  attempts  to 
observe  the  restricted  scaling  with  an  error  in  pc.  Vertical  cuts  of  die  surface  by  planes 
<T=const  give  the  results  of  parameter  variation  in  a  noise-limited  experimental 


system. 


Obviously,  the  experimental  difficulty  associated  with  the  problems 
mentioned  here  will  depend  on  the  particular  system.  Problems  will  be  minimized  if  the 
dynamic  range  of  added  noise  is  maximized  and  well  above  intrinsic  noise  levels  (for 
example,  in  the  numerical  experiments,  noise  was  present  in  the  form  of  roundoff  error, 
but  its  strength  was  separated  by  many  orders  of  magnitude  from  that  of  the  added  noise). 
This  requirement  may  conflict  with  the  avoidance  of  potential  problems  like  those 
discussed  in  Chapter  V  B,  where  the  global  properties  of  the  system  place  an  upper  bound 
on  the  noise  strength  that  can  be  meaningfully  applied. 
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VIH.  EXPERIMENTAL  CONFIRMATION  OF  SCALING  LAWS  (1)  AND  (2) 

In  this  chapter,  we  present  the  results  of  physical  experiments  that  confirm 
the  applicability  of  scaling  laws  (1)  and  (2)  to  physical  systems.  We  have  observed  crisis- 
induced  intermittency  in  a  variable-noise,  nonlinear  mechanical  system,  and  have  for  the 
first  time  confirmed  the  theoretically-predicted  key  role  played  by  an  unstable  orbit  in 
mediating  a  deterministic  crisis.  We  identify  the  unstable  periodic  orbit  mediating  the 
crisis  and  compare  the  critical  exponent  predicted  on  the  basis  of  the  orbit’s  eigenvalues 
with  the  scaling  behavior  of  the  time  between  bursts,  x.  We  find  good  agreement  between 
the  predicted  and  measured  exponents. 

We  have  also  observed  noise-induced  intermittent  bursts  in  the  nonlinear 
mechanical  system  for  p<pc,  and  have  confirmed  experimentally  the  applicability  of 
scaling  law  (2)  to  a  physical  system.  We  calculate  the  critical  exponent  y  for  the  p>Pc 
system,  using  three  operationally  independent  means  that  produce  results  in  good 
agreement,  and  show  that  the  characteristic  time  for  the  p<pc  system  is  consistent  with  a 
single  function  g(»)  over  a  wide  range  of  O' and  p-pc,  when  y  is  used  to  scale  the  results 
as  in  Eq.  (2). 

In  Section  A,  we  discuss  the  experimental  setup  and  the  measurements 
made  on  the  nonlinear  mechanical  system.  In  Section  B,  we  discuss  the  results  obtained 
by  operating  the  system  with  the  variable  noise  turned  off,  which  are  relevant  to  the 
theory  of  deterministic  crises  leading  to  Eq.  (1).  In  Section  C,  we  discuss  the  results 
obtained  by  adding  noise  to  the  system. 

A.  Experimental  Setup  and  Procedures 

The  experimental  system  was  a  gravitationally  buckled,  amorphous, 
magnetoelastic  ribbon16  (transversely  annealed  Fe8iB13.5Si3.5C2;  3mm  x  65mm  x 
25/im),  clamped  vertically  at  its  base,  and  driven  parametrically  by  a  time-varying 
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component  to  an  applied  magnetic  field  (Fig.  21).  The  time-varying  component  had  both 
a  sinusoidal  and  a  random  contribution.  The  random  component  was  given  by 

oo 

£(*)=  ^ak{e(t-kSt)-d[t-(k+l)dt]}, 

k=-oo 

where  6(*)  is  the  unit  step  function,  8t  is  an  update  interval  much  smaller  than  both  the 

mechanical  response  time  of  the  ribbon  and  the  forcing  period,  and  the  a/c  are 

2 

independent,  zero-mean  Gaussian  random  variables  with  variance  oa.  This  form 

approximated  Gaussian  white  noise.  The  ribbon  is  made  from  a  new  class  of  amorphous 
materials  that  exhibit  very  large  reversible  changes  in  their  Young’s  modulus  E  with  the 
application  of  a  small  magnetic  field  (inset,  Fig.  21).  The  oscillating  magnetic  field 
changes  the  stiffness  of  the  ribbon,  which  therefore  buckles  to  a  correspondingly  greater 
or  lesser  degree.  It  should  be  noted  that  although  the  variable  noise  appears  additively  in 
the  magnetic  field,  the  ribbon  experiences  parametric  noise;  thus,  the  applicability  of  the 
theory  leading  to  Eq.  (2)  is  not  immediately  guaranteed.  The  degree  to  which  the  ribbon 
is  buckled  is  measured  by  an  MTI  Fotonic  sensor  near  the  base  of  the  ribbon.  Additional 
details  of  the  apparatus  are  discussed  in  Ref.  17.  The  sensor  output,  monotonically  related 
to  the  ribbon  curvature,  was  used  as  the  single  measured  dynamical  variable  in  the 
experiments.  The  rest  of  the  system’s  phase  portrait  was  reconstructed  by  delay- 
coordinate  embedding.18  The  sensor  output  voltage  time  series  was  Vj  =  V(ti)  ( t{=iAt,  i=  1, 

2,  3...;  At=l/f,  the  forcing  period).  This  choice  of  At  was  made  to  obtain  a  stroboscopic 
Poincare  section  in  a  (d+ 1  )-dimensional  continuous- time  phase  space.  For  sufficiently 
large  d,  this  choice  of  At  induces  a  discrete  d-dimensional  map  xz+i  =  F(x/ )  for  the 
system,  where  x,-  =  (V) ,  V*+i,  Vj+ 2,. . .  Vi+d)-  It  should  be  noted  that  a  complete  model  of 
the  ribbon’s  behavior  requires  a  partial  differential  (i.e.,  infinite-dimensional)  equation.  It 
is  therefore  noteworthy  that  essential  features  of  the  ribbon’s  dynamics  can  be  captured  in 
a  low-dimensional  mapping,  as  has  been  shown  previously.17*19 
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FIG.  21.  Experimental  setup.  Inset:  ratio  of  Young’s  modulus  E  of  the  ribbon  to  zero- 
field  modulus  Eq,  vs  applied  magnetic  field. 
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We  observed  that  for  an  imposed  magnetic  field  H  =  Hdc  +  H acsin(2  jrft), 
(where  Hdc  =  0.82  Oe  and  Hac  =  2.05  Oe),  as /decreases  through  about  0.97  Hz,  the 
“deterministic”  (<f^0)  system  undergoes  a  crisis.  [It  is  important  to  note  that  the  crisis 
occurs  as /decreases  through  its  critical  value,  which  is  the  opposite  of  the  convention  for 
the  generic  crisis  parameter  p  used  previously  in  this  dissertation  and  in  many  of  the 
references.  To  minimize  confusion,  we  use  /  explicitly,  rather  than  p  in  this  chapter.] 
Before  the  crisis,  the  system  state  moves  chaotically  on  a  strange  attractor  [Fig.  22(a)]. 
After  the  crisis  [Fig.  22(b)],  the  attractor  is  much  larger.  The  dynamics  during  the  bursts 
is  more  complicated  than  that  between  bursts.  However,  the  part  of  the  attractor 
corresponding  to  the  pre-crisis  attractor  (the  core)  is  adequately  represented  in  a  two- 
dimensional  phase  space,  and  so  we  expect  the  theory  leading  to  Eqs.  (4)  and  (5)  to  apply. 
Turning  on  the  noise  term  (§*0)  produces  orbits  much  like  those  of  the  post-crisis  system, 
even  for  f>fc  [Fig.  22(c)]. 
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V,  Vi 


FIG.  22.  Delay  coordinate  embedding  of  time  series  taken  before  (a)  and  after  (b)  the 
crisis.  The  time  delay  At  is  1  If,  the  stroboscopic  sampling  period.  The  core  attractor  of 
(a)  is  enlarged  by  burst  dynamics  in  (b).  In  (c),  the  parameter  is  held  at  the  same  value 
as  in  (a),  but  the  variable  noise  is  on.  The  noise  produces  burst  dynamics  similar  to 
those  seen  in  (b).  The  small  crosses  in  (a)  indicate  the  elements  of  the  period-3 
unstable  orbit  that  mediates  this  crisis. 
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B.  Analysis  of  Deterministic  Crisis 

The  scaling  exponent  y  of  the  deterministic  crisis  was  calculated  in  three 
operationally  independent  ways.  First,  y  was  determined  by  measuring  the  characteristic 

A 

time  between  bursts  r  as  a  function  of  f-fc,  and  fitting  the  results  to  Eq.  (1).  Estimates  x 
of  t  were  obtained  from  3-hour  time  series  recorded  for  closely  spaced  values  off  near 
the  crisis  (T  is  measured  in  forcing  periods  of  the  sinusoidal  magnetic  field  component). 
The  theoretical  scaling  law  (1)  for  x  assumes  that  the  system  is  deterministic;  for  such  a 
system,  the  characteristic  time  goes  to  infinity  at  the  crisis.  Real  systems,  which  are 
always  accompanied  by  noise,  show  a  finite  characteristic  time  at  f=fc  as  is  evident  from 
Eq.  (2).  The  presence  of  residual  system  noise,  even  with  the  a/  set  to  zero,  made 
determining  the  value  of  fc  more  complicated  than  in  an  ideal  system,  and  somewhat 
limited  the  accuracy  with  which  it  could  be  determined.  The  value  of  fc  for  this  crisis  was 
determined  by  minimizing  the  degree  to  which  a  plot  of  x  vs  log (fc~f)  deviated  from  a 
straight  line,  using  data  taken  not  too  near  fc  [where  noise  more  strongly  affects  the 
applicability  of  Eq.  (1)].  Given  fc,  y  was  determined  from  the  best  nonlinear  least-squares 
fit  of  the  functional  form  x  =  k(fc~f)~V  to  the  data.  The  data  and  corresponding  fit  are 
shown  in  Fig.  23.  The  statistical  estimate  of  the  critical  exponent  using  this  method  is  y= 
1.12±0.02. 

A  second  determination  of  y  was  made  by  considering  the  accumulation 
of  probability  measure  at  the  edge  of  the  attractor.  The  theory  leading  to  scaling  law  (1) 
assumes5  that  the  measure  fi  of  the  attractor  within  e  of  its  edge  accumulates  as  ~  eX 
Using  a  trajectory  of  several  hours  duration  recorded  for /just  above  fc,  the  iterates  were 
sorted  in  magnitude,  and  an  estimate  of  f±(£)  was  obtained  as  a  cumulative  count  of  the 
number  of  iterates  within  e  of  the  maximum  iterate.  The  slope  of  the  straight-line  portion 
of  a  plot  of  ln[/i(e)]  vs  ln(e)  produced  an  estimate  of  y=  1.08±0.05,  which  is  statistically 
consistent  with  the  above  estimate. 
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A  final  determination  of  y  resulted  from  the  observation  that  the  crisis 
shown  in  Fig.  22  is  a  homoclinic  tangency  crisis  of  a  two-dimensional  map.  Careful 
observation  of  the  system  state  just  before  a  burst  indicated  that  system  trajectories 
closely  approached  a  period-3  orbit  [the  elements  of  which  are  indicated  in  Fig.  22(a)] 
before  leaving  the  core  attractor.  Figure  24  shows  V/+ 3  vs  Vj  for  the  phase-space  region  of 
the  core  attractor  after  crisis;  several  iterates  passing  near  the  saddle  are  indicated.  Any 
period-three  orbit  appears  as  fixed  points  of  the  three-times-iterated  mapping  F(^)  and, 
because  we  use  a  delay-coordinate  embedding,  these  must  lie  along  the  45°  line  (i.e., 
Vi+ 3  =  Vj).  The  attractor  intersects  the  45°  line  in  four  places.  One  of  these  intersections  is 
a  period- 1  point;  the  other  three  are  elements  of  the  same  period-3  saddle.  For  f>fc  (before 
the  crisis),  the  attractor  only  intersects  the  45°  line  at  the  period- 1  point.  This  indicates 
the  absence  of  any  period-3  orbits  on  the  attractor  before  crisis.  Crises  in  two-dimensional 
maps  are  limited  to  two  varieties,  depending  on  whether  the  mediating  periodic  orbit  is  on 
the  attractor  before  (heteroclinic),  or  only  just  at  and  after  (homoclinic)  the  crisis.^  The 
absence  of  a  period-3  orbit  on  the  attractor  before  the  crisis  means  that  the  crisis  is  of  the 
homoclinic  type.  The  theoretically  predicted  critical  exponent  for  this  crisis  should 
therefore  be  given  in  terms  of  the  eigenvalues  of  the  period-3  orbit  by  Eq.  (5). 

Estimating  these  eigenvalues  from  the  experimental  data  is  complicated  by 
the  facts  that:  (i)  the  map  F  is  clearly  very  dissipative,  (it)  the  system  is  contaminated  by 
noise,  and  (iii)  the  linear  neighborhood  of  the  period-3  point  is  relatively  small. 
Therefore,  following  an  orbit  as  it  passes  near  the  period-3  point,  as  suggested  in  Ref.  5, 
is  not  a  good  way  to  estimate  the  eigenvalues.  Instead,  we  use  pairs  of  iterates  to  estimate 
the  Jacobian  matrix  of  the  map  around  the  period-3  point,  and  from  this  estimate  obtain 
the  linearized  eigenvalues  of  the  periodic  point  [see  Appendix  C].  The  estimate  of  the 
critical  exponent  based  on  these  eigenvalues  of  the  period-3  saddle,  y—  1.07±0.15,  clearly 
is  statistically  consistent  with  the  exponents  determined  independently  above.  Such 
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consistency  strongly  supports  the  key  role  played  by  unstable  periodic  points  in 
mediating  crisis  phenomena  in  nonlinear  dynamical  systems. 


C.  Analysis  of  Noise-induced  Crisis 

Although  the  analysis  in  the  previous  section  was  carried  out  as  though  the 

system  operated  with  0  were  deterministic,  even  then  there  is  some  residual  noise  in 

the  system,  resulting  from  mechanical  vibrations  and  imperfections  in  the  operational 

amplifier.  This  residual  system  noise  must  be  taken  into  account  when  applying  Eq.  (2), 

in  which  the  strength  <7  of  the  total  noise  experienced  by  the  system  appears.  We  allowed 

for  this  fact  as  follows.  When  the  system  was  operated  at  a  different  set  of  parameters 

(Hdc  =  0.22  Oe,  Hoc  =  2.05  Oe,/=0.95  Hz)  it  was  noted  that  the  attractor  was  not  chaotic, 

but  rather  was  a  stable  periodic  orbit.  The  sequence  of  iterates  Vi  had  a  Gaussian 
distribution,  due  to  the  system  noise,  with  a  standard  deviation  <\  of  about  5mV.  As  the 

variable  noise  was  increased,  the  distribution  remained  gaussian,  but  broadened.  As 
shown  in  Fig.  25  the  effect  of  the  variable  noise  on  the  system  state  was  well  described 
by  an  additive  variance  model,  consistent  with  the  residual  system  noise  and  the  variable 

noise  being  independent  Gaussian  processes.  In  testing  the  applicability  of  Eq.  (2),  we 
used  an  effective  noise  strength  <7=  [Oq  +  cra]  , where  <tq,  the  residual  system  noise, 

was  determined  from  the  fit  in  Fig.  25. 

Figure  26(a)  shows  the  characteristic  times  T  estimated  from  data  runs 
taken  at  a  wide  variety  of  noise  and  (pre-crisis,  f>fc)  parameter  values.  In  Fig.  26(b),  the 
scaling  indicated  by  Eq.  (2)  has  been  applied,  using  the  critical  exponent  determined  for 
f<fc  and  a/c= 0.  The  widely  dispersed  data  of  Fig.  26(a)  now  collapse  to  a  single  curve, 
which  gives  the  graph  of  the  function  g(»)  in  Eq.  (2)  for  this  system.  This  consistency 
strongly  supports  the  theory  underlying  scaling  law  (2). 
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FIG.  25.  Measured  standard  deviation  of  the  iterate  sequence  Vi,  as  a  function  of  the 
strength  of  the  added  noise.  Here,  the  system’s  attractor  is  a  stable  periodic  orbit.  The 
°  2  2  1/2 
curve  shows  the  best  fit  to  the  sum-in-quadrature  model  av  =  Ic[<Jq  +  <7a\ 
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FIG.  26.  Characteristic  times  for  noise-induced  bursting,  (a)  shows  range  of  raw  data. 
Scaling  of  Eq.  (2)  is  applied  in  (b),  using  deterministic  critical  exponent  7.  Data 
collapse  onto  a  single  curve,  giving  graph  of  function  g(*)  in  Eq.  (2). 
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IX.  CONCLUSIONS 


In  this  dissertation  we  have  introduced  a  general  scaling  law  for  the 
characteristic  time  T  of  noisy  systems  operating  at  parameter  values  near  those 
corresponding  to  a  deterministic  crisis.  The  scaling  law  (2),  applicable  to  systems  whose 
deterministic  crises  are  governed  by  the  scaling  law  (1),  shows  that,  although  the  strength 
of  noise  in  the  system  is  fundamental  in  determining  the  characteristic  time  T,  the  noise- 
induced  crisis  inherits  its  governing  exponent  y  from  the  deterministic  system.  For  the 
case  of  the  noisy  logistic  map,  our  results  reduce  to  a  previously  known  special  case.9 

Numerical  examples  suggest  that  the  noise  scaling  law  can  be  expected  to 
apply  to  crises  resulting  in  attractor  destruction,  attractor  enlargement,  and  attractor 
merging,  for  both  maps  and  flows.  However,  some  complications  to  the  observability  of 
the  scaling  law  have  also  been  identified. 

A  physical  experiment  further  supports  the  theory  underlying  scaling  laws 
(1)  and  (2).  The  fact  that  the  noise  in  the  experiment  appeared  parametrically  shows  that 
the  applicability  of  Eq.  (2)  is  even  wider  than  suggested  by  the  numerical  results. 

The  ubiquity  of  noise  in  physical  systems  means  that  the  experimental 
investigation  of  systems  near  crisis  should  be  done  with  congnizance  of  scaling  law  (2). 
Conversely,  the  inheritance  of  the  critical  exponent  from  the  deterministic  crisis  by  the 
noise-induced  crisis  means  that  the  intentional  addition  of  noise  may  be  an  important 
experimental  probe  for  investigating  the  deterministic  dynamics  underlying  a  noise- 
limited  system. 
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APPENDIX  A:  DERIVATION  OF  EQS.  (12)  AND  (13) 

In  this  appendix  we  derive  the  first-order  correction  terms  to  the  overall 

rate  of  escape  for  orbits  of  the  noisy  logistic  map,  Eqs.  (12)  and  (13).  The  left-hand 
escape  rate  L  is  derived  by  first  considering  the  dynamical  density  or^ts  near 

the  right  edge  of  the  attractor,  after  noise  has  been  added,  but  before  the  deterministic  part 
of  the  map  has  again  been  applied: 

oo 

Px+Z(w)  =  J  p(x)  P^w-x)  dx 

— OO 

p! 4 

'  1  1 
crJlK 


exp[-{w-x)2l2(?\  dx. 


Making  the  substitutions  s  =  w-p/4  and  u  =  p/4-x,  we  have 


/• 

,.-1/2  „„„ 

/  2  A 

VT  _S_ 

u  1  exp 

«/ 
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and  using  an  integral  definition  of  the  parabolic  cylinder  function  [see  Ref.  20,  Eq.  3.462 
(1)] ,  we  get 

Px+fi')  =  D-l/2(s/o). 


The  integral  of  this  result  over  [1,°°)  yields  the  right-hand  escape  rate  R. 
By  applying  the  deterministic  part  of  the  map  (12),  one  obtains  the  density  of  orbits  near 
the  left  side  of  the  attractor,  before  the  noise  kick: 

Px+Z(x)  =  3/2  \cr^2eXP^  (1-P/4-x/4)2/4°21 

x£>-1/2[(  l-p/4~x/4)lo] . 

The  density  of  orbits  after  the  noise  is  added  is  obtained  by  convolving  the 
noise  density  with  the  part  of  P^+Z W  to  t^e  right  of  x=0.  The  integral  of  the  resulting 

convolution  gives  the  left-hand  escape  rate  L.  Because  the  convolution  does  not  exist  in 
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closed  form,  we  use  a  polynomial  approximation  to  P %+z(x)  in  the  neighborhood  of  x=0 
(we  use  an  asymptotic  expansion): 


'-max* 


PX+2^  ~  X  lc\ ' 

k= 0 

Now  the  convolution  integral  can  be  evaluated: 

9  T  kmax~(k)  xk 

L  -  jdy  jdx  tpxlz(^pz(y-x^ 
0  k= 0 


or,  substituting  u=y-x,  we  have 

9  Y  *5 max ~  ,,,  (y-u)k 

=  jdy  jdu  ^  P  x+zfQ'id  Pz9^' 
k= 0 


Reversing  the  order  of  the  sum  over  k  and  the  integrations,  expanding  the  binomial, 
integrating  by  parts,  and  letting  w=  -y,  we  have 
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Using  the  definition  of  the  gamma  function,  we  have 


?f+z(0) 

kjnax  Jc+Xyik-Y)!! 
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>< 

7=0  7 


The  sum  over  j  can  be  simplified,  yielding 

\k 


J.  1 


PxlzfO) 
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l=  X  -nr- (-»*-  —i/2  n^i)(-i),‘faT 

£=0  ^ 
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Finally,  we  use  the  following  recursion  relation  for  the  parabolic  cylinder  functions, 

J^t)=fD*00-Afc+ito, 

(see  Ref.  12,  Eq.  19.6.2)  to  find  the  derivative 

?xlz« d-P/^MO2] 


yielding 


x  ^P-l/2[(  1  — p/4— jt/4)/<r] , 

kmax 


L^Ar^\  EoST&D‘-'«[(1-p/4VCTl  ■ 

k= 0 

Using  the  fact  that  for  the  logistic  map  pc= 4,  we  obtain  Eq.  (12). 

The  steps  in  deriving  R ,  the  failure  rate  for  right-hand  escapes,  are  exactly 
the  same,  except  the  limits  of  integration  are  permuted: 

OO  0  k  I 

_  ,  ,  Kmax~(l/s  xk 

R  =  U  J<fe  %Pf^0)^P^y-x). 

0  -oo  /:=0 


The  result  is  the  same  as  for  L,  except  that  the  series  is  an  alternating  one,  yielding  Eq. 
(13). 
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APPENDIX  B:  INTEGRATION  OF 
STOCHASTIC  DIFFERENTIAL  EQUATIONS 


Obtaining  realizations  of  stochastic  differential  systems  such  as  Eq.  (17)  is 
a  less  well-developed  branch  of  numerical  analysis  than  the  numerical  solution  of 
ordinary  differential  equations.  Various  ad  hoc  schemes  have  appeared  in  the  literature 
(cf.  Refs.  8  and  9)  that  undoubtedly  provide  qualitatively  correct  results,  but  which  may 
not  be  accurate  enough  when  numerical  results  are  to  be  compared  with  experimental 
results,  or,  as  discussed  below,  with  numerical  results  for  the  corresponding  deterministic 
system.  For  the  results  presented  in  Chapter  VI,  we  used  a  one-step  colocation  integration 
scheme  closely  related  to  that  reported  in  Ref.  21.  The  idea  of  the  method  is  to  use  a 
single  pseudorandom  vector  to  account  for  the  cumulative  effect  of  the  stochastic  process 
£(r)  over  the  whole  integration  time  step,  respecting  the  influence  of  the  dynamics  on  the 
noise-perturbed  orbit  throughout  the  time  step. 

The  method  we  used  has  an  associated  single-step  error  which  is  o(At 5/2) 
in  the  integration  time  step  At.  This  is  much  larger  than  the  single-step  error  associated 
with  most  numerical  methods  widely  used  to  solve  deterministic  ODEs.  The 
characteristic  behavior  of  random  walks  (displacement- At^ty  means  that,  for  the  n*  level 
of  successive  approximation  in  the  numerical  method,  the  error  terms  are  o(Atn '*). 
Unfortunately,  the  complexity  of  the  the  stochastic  integrals  arising  in  the  numerical 
method  becomes  prohibitive  very  quickly,  and  errors  of  order  o(At^)  are  all  that  have 
ever  been  used  in  practice.  Further,  the  use  of  adaptive  step-size  control  in  a  stochastic 
integrator  is  very  difficult  (and  no  such  algorithms  have  been  published). 

The  problem  caused  by  the  large  error  is  manifested  when  simulated 
stochastic  results  are  to  be  compared  to  numerical  results  from  the  corresponding 
deterministic  system.  A  numerical  integration  scheme  is  really  a  discrete  map,  where  one 
of  the  parameters  is  the  integration  time  step.  Thus  even  a  deterministic  system. 
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integrated  by  two  different  schemes,  or  by  the  same  scheme  using  two  different  time 
steps,  may  well  produce  different  results.  In  particular,  the  parameter  values  at  which 
crises  occur  are  likely  to  vary,  depending  on  the  integration  scheme.  In  order  to  avoid  any 
such  problems  with  our  comparisons  of  the  deterministic  and  noisy  Duffing  oscillators, 
we  integrated  the  deterministic  system  using  the  exact  same  algorithm  and  time  step  as 
for  the  noisy  system,  but  omitting  the  noise  terms.  The  resulting  second-order  Euler 
algorithm  yielded  numerical  results  that  differed  slightly  from  those  previously  reported^ 
for  the  deterministic  version  of  Eq.  (17). 
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APPENDIX  C:  ESTIMATION  OF  EIGENVALUES  OF  UNSTABLE  SADDLES 

FROM  EXPERIMENTAL  DATA 

Our  general  eigenvalue  estimation  procedure  is  outlined  as  follows.  To 

obtain  the  eigenvalues  of  a  period- «  orbit,  we  approximate  the  map  F(")  (the  n-times- 

iterated  map  F)  by  a  Taylor  series 

xi+n  +  dxi+n  =  F  («)(x/)  +  DF(«)IX.  •  dxi  +  OQSxfi), 

and  let  x;  be  the  period-n  point,  which  we  denote  x*.  Because  F(n)(x*)=x*,  we  see  that 
the  Jacobian  DF(")IX*  takes  tangent  vectors  into  tangent  vectors  near  the  saddle.  We  can 

use  the  difference  between  a  pair  of  iterates  on  the  same  orbit  to  approximate  a  tangent 

vector,  and  watch  its  propagation  near  the  period-n  point: 

x,+n  -  Xj+n  «  DF(«)IX*  •  (x;  -  xj),  i  *  j. 

The  key  advantage  of  this  method  over  that  of  Ref.  5  (1987)  is  that  here  we  can  take 
advantage  of  widely  (temporally)  separated  recurrent  iterates  (i.e.,  i  very  different  from./). 
Furthermore  we  do  not  require  many  sequential  iterates  which  are  all  within  the  linear 
neighborhood  of  the  saddle  and  whose  differences  are  also  significantly  above  the  system 
noise  level. 

Because  we  use  a  delay-coordinate  embedding,  the  surface  of  section  map 

for  x=(y,z )  is  of  the  form 

iyi+n  =  h(yi,Zi) 

{Zi+n=yi 

and  we  therefore  know  that  the  form  of  the  Jacobian  is  particularly  simple: 


where  r=dh/dy\^  z*^  and  s=dh/dz\^^  z*y 

Calculating  the  eigenvalues  of  DF(n)lx*  explicitly  and  using  Eq.  (5)  and 
the  fact  that  /Ji/fe  =  det^  DFWlx*  we  have  y=  log[r/2  -  \(r  2+  As)l^\!\_2  loglsl)].  Thus, 
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estimation  of  the  eigenvalues  of  the  period- n  orbit  is  reduced  to  the  statistical  estimation 
of  the  quantities  r  and  s,  from  the  linear  multiple  regression  model 
yi+n  ~  yj+n  =r(yt-  yj)  +  s  (z;  -  zj) , 

where  a  set  of  (i  ^  j)  pairs  is  chosen  from  the  vicinity  of  the  saddle.  Constraints  can  be 
placed  on  the  regression,  because  the  sign  of  s  will  be  known  from  the  type  of  period-^ 
saddle  involved  (an  ordinary  saddle  implies  s<0;  a  flip  saddle  implies  s>0).  Furthermore, 
the  (i  &  j )  pairs  can  be  chosen  in  such  a  way  as  to  maximize  the  statistical  confidence  of 
the  estimates,  i.e.,  chosen  such  that  the  contributions  of  the  r  and  s  terms  in  the  regression 
are  about  equal.  [Use  of  the  restricted  form  of  the  Jacobian  and  a  careful  choice  of  data 
for  the  regression  analysis  is  essential  in  a  noisy  system.  Least-squares  estimation  of  the 
general  Jacobian  with  a  comprehensive  data  set  was  done  successfully  with  a  lower-noise 
system  in  Ref.  22,  but  for  our  system,  this  results  in  a  very  poor  estimate  of  the  stable 
eigenvalue,  and  consequently,  a  poor  estimate  of  y.]  Finally,  propagation  of  errors  using 
the  y(r,s)  formula  allows  the  statistical  uncertainties  in  r  and  s  to  be  translated  into 
uncertainty  in  y  A  nearly  identical  analysis  pertains  to  the  heteroclinic  case. 

For  the  particular  crisis  of  interest  in  Chapter  VIII,  we  take  the  point  x*  to 
be  as  shown  in  Fig.  4.  Observation  of  iterates  passing  near  x*  shows  that  the  period-3 
saddle  is  an  ordinary  one  (s<0),  that  the  stable  direction  is  nearly  vertical ,  and  that  the 
unstable  direction  has  a  slope  of  about  0.1.  The  second  row  of  DF  indicates  that  the 
eigenvectors  must  be  of  the  form  (/J,l).  We  therefore  know  that  p2«l,  while  /?i~10, 
indicating  that  lsl«lrl,  i.e.,  DF  is  ill-conditioned.  This  information  allows  the 
determination  that  the  best  regression  estimates  of  r  and  s  will  be  obtained  for  pairs  of 
iterates  whose  vector  difference  is  nearly  vertical.  Using  150  such  pairs  from  a  single 
orbit,  we  estimated  r,  s  and  their  statistical  uncertainties,  and  conclude  that  y-  1.07±0.15. 
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